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ABSTRACT 


The dihedral angle distribution is a useful parameter in the quantification of 
polycrystalline microstructures. However, this microstructural parameter has received not 
much attention because the methodologies associated with the accurate measurements of 
dihedral angles are generally tedious and time consuming. In this study an automated 
methodology based on image processing and analysis, was developed to evaluate dihedral 
angles in microstructures in a two-dimensional plane section of polycrystals It is shown 
that the technique is extremely fast and accurate and therefore can be applied to obtain a 
significantly large statistical sample of dihedral angles. 

The study was expanded further to transform the plane dihedral angle distribution to true 
dihedral angle distributions using a 3-dimensional probe termed as “disector”. The true 
dihedral angles were further related to grain boundary energies to characterize 
polycrystals. This probe introduces precision as well as efficiency over the classical 
single as well as serial sectioning technique. This technique gives more information than 
a single 2-dimensional probe and is far less tedious and less complicated than the serial 
sectioning technique 

The application of the above methods for calculating plane and the true dihedral angles is 
demonstrated by evaluating triple junction geometry of polycrystalline 316L austenitic 
stainless steel. The validity of the measurements was established by a comparison of its 
result with careful manual measurement of results. 
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CHAPTER 1 


INTRODUCTION AND OBJECTIVES 


Interfaces are a fundamental constituent of the polycrystalline microstructure. The 
geometry of these interfaces is associated with their energy state, which is strongly 
linked with mechanical properties of materials (Singh 1998). In a polycrystalline 
space filling microstructure mating of three grains form a triple edge. The dihedral 
angles are the angles between any two surfaces (or interfaces) that meet at the triple 
edge. These angles have been termed as true dihedral angles (Murr 1975). While , in 
a polycrystalline microstructure observed on a two-dimensional section, the angles 
between grain boundary segments at a triple point are termed as plane dihedral 
angles. The geometry of triple junctions, described in terms of the true dihedral 
angles between grain boundaries, is related to the energy state of grain boundaries 
and consequently to the deformation behavior of polycrystals [12,14-16]. Figure 1 
shows the description of triple point geometry by the angles between the grain 
boundary segments on a polycrystalline micrograph of a plane section. 



2 



Figure 1: An optical microstructure of a polycrystallme material, pi, [32 
and Ps are the plane dihedral angles at a triple point. 


Even though dihedral angles measure has been considered to be an important 
parameter for the characterization of polycrystals, it has been used to characterize 
polycrystalline microstructure only to a limited extent . This is attributed to the fact 
that the procedures generally followed are tedious and the results may not reflect 
good statistical estimate of the distribution of plane dihedral angles. Therefore the 
main objective of this study is to develop an automated technique based on digital 
image processing and analysis to obtain a large and an accurate statistical estimate of 
the distribution of plane dihedral angles. For testing the validity of the technique it 
has been applied to samples of 316L austenitic stainless steel. 


For investigation of triple junction geometry numerous methodologies are available 
among which quantitative stereology promises to be very efficient. It provides a 
means to describe the geometry of three-dimensional microstructure from 
observations made on plane two-dimensional sections. Classical stereological 
methods typically involve measurements made on a single plane section. The 
disadvantages of using a single section are well known (Gundersen 1986) and can be 
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completely overcome by using the technique of serial sectioning. However, the serial 
sectioning technique is very cumbersome and time consuming. Some of the 
disadvantages of both the single section and serial sectioning techniques have been 
overcome by the three-dimensional sampling probe termed as the “disector” (Sterio, 
1984, Gundersen, 1988). The disector technique, which uses two parallel plane 
sections separated by known distance as shown in figure 2, has been used for an 
unbiased determination of numerical density, Ny (Gundersen 1986) 



Figure 2: Illustration of two parallel section of a polycrystal and true 
dihedral angles. 

(Murr 1975) 



The other objective is to directly evaluate true dihedral angles. It is shown that a 
statistically significant number of measurements of true dihedral angles can be 
performed using automated digital image processing. Further, these measurements 
have been used to obtain a distribution of relative grain boundary energies. 

Therefore to summarize, main objective of this study is to develop a combination of 
an automated technique based on digital image processing/analysis and the three- 
dimensional probe, disector, to obtain a large and an accurate statistical estimate of 
distribution of plane dihedral angles from two sections of the given sample and to 
transform it to true dihedral angle distributions. Further, to show that the true dihedral 
angles can be used for evaluating relative grain boundary energies. 
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CHAPTER 2 

MICROSTRUCTURAL CHARACTERIZATION OF 
POLYCRYSTALS: AN OVERVIEW 


2,1 Characterization of Polycrystals 


To characterize a polycrystalline microstructure various parameters, like grain size 
distributions, grain volume, surface area measurements, curvatures etc, have been 
extensively examined by studying the two-dimensional planar sections using quantitative 
stereology. 


2.1.1 Classical Stereology 

Stereology has emerged as an effective tool to characterize a microstructure. Quantitative 
stereology is generalized body of methods, applicable to any material, for quantitative 
characterization of three-dimensional microstructures by means of two-dimensional 
sections taken at random from the solid material (Underwood, 1970). Much of the earlier 
works since has concentrated on interpretation of two-dimensional measurements 
obtained from a planar section and relating it to three-dimensional parameters defining 
the structure (Wiebel, 1980). 
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The methods involved in stereology allow us to infer parameters of the structure from 
incomplete information gathered from sections of a dimension lower than that of the 
structure This has two consequences; firstly certain parameters of structure (connectivity 
between parts) are not accessible for stereological evaluation and secondly the section 
being a sample of a particular type, stereological methods becomes statistical in nature 
and measurements obtained on these seimples yield estimates, which are averages and are 
affected by some sampling error. This limitation requires that the unit on which the 
parameter is defined occur in large number within the structure. Thus it is easy to 
estimate the total surface, average diameter and the size distribution of a population but it 
is impossible to estimate the surface or the diameter of a single object by stereological 
methods, unless we are able to cut a large number of random sections from this one and 
only object. A section or probe is a statistical sample of the structure, which is 
independent of the geometry of the structure. The section image that we observe must be 
a chance intersection between structures and probe i.e. random sampling. In 1945 
Saltykov showed that the boundary length of profiles on sections was proportional to the 
surface area of the embedded objects and this surface can be estimated from the number 
of intersections the boundary form with a random test line. It was also observed that the 
principle establishing the probabilistic relationship between the area of a surface and 
length of a line, expressed by the number of intersections formed when they interact in 
random fashion as shown in figure 2.1. It can be seen that Nl= 5 (number of line 
segments formed) and Pl= 5 (number of points formed by intersection of the test line and 
the grain boundaries), then the surface area of the grain boundary in a unit volume is 
given by: 


Sv = 2Pl = 2Nl (2.1) 

This principle was also extended to estimate the length of linear features of the 
structures like grain boundary edges. The reference system is now a set of planes instead 
of test lines and the observer counts the number of times the lineal features transect the 
plane of defined area. 
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The fundamentals of stereological symbols are essentially noted as a double symbol of 
two capital letters. The first capital letter defines the component parameter and the second 
written as subscript, represents the parameter of the reference space. To illustrate the 
symbol; the numerical density of some particulate component (b) in the containing space 
(c) is written as 

Nv =N(b)/V(c) 

And the areal density of profiles of (b) in the profile of structure (c) is represented as 
AA = A(b)/A(c) 



polycrystalline microstructure, intersecting 
grain boundaries at five locations 
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DISTRIBUTIONS 

Stereo logical principles were also developed to explicitly consider the individual object 
or element of structure while determining the number and size of particles. Sectional 
profiles are mostly smaller than the object they have been derived. Stereological methods 
are now available by which we can transform measurements of profile size into 
dimensions describing the size of the objects. These problems were first tackled by 
Wicksell (1925) to solve a corpuscular problem in anatomy. He derived relationship 
between the size distribution of spherical and ellipsoidal particles with that of their 
profiles on the section. Saltykov (1967) developed a method for transforming distribution 
of two-dimensional measurement of sections to three-dimensional distribution of spheres. 
The problem was complicated by the fact that spheres with different diameters j, j-i-1 , jH-2 
etc will contribute sections of the same size ‘i’, so source of each section were 
determined. To illustrate the transformations for size distributions a poly dispersed system 
of spheres that had been separated into five class intervals. The number of spheres per 
unit volume in each class interval is 

(Nv)i, (Nv)2, (Nv)3. (Nv) 4 and (Nv)5 
and their corresponding radii are 


ri,r 2 , rs, r4 and rs 


the total number of circular sections of a given size per unit area, regardless of origin, is 
(Na)i, (Na) 2 , (Na)3, (Na) 4. and (Na) 5 
where the section radii vary fi:om 


0 to ri, ri to r 2 , r 2 to rs, rs to r 4 and xa to rs 
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Figure 2.2 illustrates the possibilities arising from a random plane intersecting particles of 
all five sizes. Spheres of the largest class interval will give sections of all radii when cut 
by a random plane also the observed number of sections of the smallest size per unit area, 
(Na)i, is sum of contributions from spheres of all sizes represented as: 


(Na)i= E (Na),j = (Na)i.i + (Na)i,2 + (Na)i,3 + (Na)i,4 + (Na)u (2.2) 


where i represents constant section size and j represents the sphere diameters. The largest 
section on the plane of polish can unambiguously belong to the largest sphere only which 
is represented as: 


(Na)5,5 = (Na)5 

(Na)4,4 = (Na)4 - (Na)4,5 

(Na)3,3 = (Na)3 - (Na)3,5 - (Na)3,4 

(Na)2,2 = (Na)2 - (Na)2,5 - (Na)2,4 - (Na)2,3 

(Na)u = (Na)i - (Na)i,5 - (Na)i,4 - (Na)i,3 - (Na)i,2 (2.3) 


(Na) 5 , is measured from the test plane , (Na) 5,5 is known 

Further in order to deduce how many sections of each size come from spheres in each 
class interval, probabilities of random planes intersecting a sphere to give sections of 
each size. If P,j is the probability of a test plane intersecting a sphere within a slice of 
height h, (as shown in figure 2.3) is given by: 
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/j/r„ 


h,_^ -h, I r„ 


(2.4) 


P, = (2.5) 

^max 

if the polydispersed system of spheres is divided into class intervals then, 
(.K),='Z(N,),,,ID^ mm-’ 

I 


'LiI^X,=P,DAP^y), (2-6) 

J 

(Na) 5, is measured from the test plane , therefore (Na) 5,5 is known from equations 2.3. P,,, 
is calculated from equation 2.5; D5 is measured from the plane of polish. (Nv)5 is thus 
evaluated from equation 2.6. Knowing (Nv)5 we proceed to calculate (Nv)4, from 
equation 2.5 probability of hitting size-5 sphere to get a size-4 section. This value of P4,5 
in equation 2.6 gives which when subtracted from the measured value (Na) 4, 

gives (Na)4,4 from equation 2.3. Using the appropriate probability for hitting a size-4 
sphere to get a size-4 section and knowing (Na) 4,4 , (Nv)4 is calculated using equation 2.6. 
These sequences of calculations are continued down to the smallest particle size. It is 
seen that each subtracted term depends upon the previously calculated values of number 
density, and that the number of the smallest particle maybe inaccurate because of 
accumulated errors. 
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of spheres to give sections of 

diameters di to ds ^ ^ ^ 

Figure 2.3: geometrical relationship in 

the intersection of a sphere by a plane 

with a slice h. 


2.1.2 Grain boundary energy as a parameter 


Another parameter to characterize polycrystals is grain boundary energy associated with 
the geometry of the junctions where three grain boundary planes meet. This parameter 
can be used to determine properties like plastic deformation of any polycrystalline 
material (Singh and Sangal, 1998). The grain boundary energy distribution has been 
investigated by Phumphery et al. (1974) and Varin et al. (1980) by studying the spreading 
and disappearance kinetics of grain boundary locations in TEM samples However the 
technique is cumbersome and prone to error in the estimation temperature and time of 
disappearance of EGDB. Murr et al. has evaluated grain boundary energy ratio by 
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directly measuring true dihedral angles using TEM images in which the three 
dimensional structures of a grain boundary planes meeting at triple edges can be observed 
as shown in figure 2.4a. In a polycrystalline microstructure, the plane dihedral angles are 
the angles at the triplepoints as indicated in figure 2 4c. Even though dihedral angles 
have been considered to be important parameters for the characterization of polycrystals, 
they have been used to characterize polycrystalline microstructure only to a limited 
extent. This is attributed to the fact that the procedures generally followed are tedious and 
the results may not reflect good statistical estimate of distribution of plane dihedral 
angles. Sangal (1998) has evolved alternate approaches to estimate distribution of true 
dihedral angles more efficiently from plane dihedral angles using rigorous stereological 
analysis. 



Figure 2.4a: illustration of geometry 
of three planes of the grains 
meeting at common junction. 

and are the set of 
dihedral angles. 



PLANE OF POLl 



/ 

/ 


POLYCRYSTAL 

STRUCTURE 

(X) 


Figure 2.4b: A polycrystalline structur 
(X) consisting of randomly oriented 
triple edges being sectioned by the 
plane of polish (Ti). 



Figure 2c: Pi, P 2 and P 3 are the plane dihedral angles at a 
triple point. 
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Estimation of true dihedral angle (TDA) distribution 

In a polycrystals the grain boundaries may have different values of energy, which may 
lead to difference in distribution of TDA (angle between the grain boundary planes 
meeting at a common edge, knovm as triple edge). TDA can be measured as mentioned 
above by TEM (transmission electron microscope) but were found to be very tedious and 
it is not possible to estimate statistical distribution of TDA. In practice the angles are 
measured as plane dihedral angles (PDA) between grain boundary segments meeting at a 
triple point. In a polycrystal, the triple edges are randomly oriented with respect to the 
plane of polish as shown in figure 2.4b. Therefore, the distribution of PDA has two 
components (1) variation in the TDA distribution and (2) statistical variation due to 
random orientation of the sectioning plane. The statistical variation tends to overshadow 
the actual variation in the TDA distribution. Hence such a distribution of PDA cannot be 
directly used for the calculation of relative energy of the grain boundaries. On the other 
hand, these surface measurements are relatively easy and it is possible to make a large 
number of measurements in a short time Therefore, in the following sections, a 
transformation procedure is developed to unfold the distribution of true dihedral angles 
from the measured distribution of plane dihedral angles. 

Relation between TDA and PDA: 

Figure 2.5 shows a randomly oriented triple edge sectioned by the plane of polish, P. the 
true dihedral angle between the grain boundary planes A and B is a and the 
corresponding grain boundary segments a and b make a plane dihedral angle of f on the 
plane of polish as shown in figure 2.5. The triple edge is parallel to the z-axis and the 
grain boundary plane A is in the XZ plane. The orientation of the plane of polish is 
described by the orientation of the plane normal (denoted by «) in terms of angle 
(f) and 6 in spherical coordinates. The unit vectors r, and Tj normal to the grain 

boundary planes A and B respectively. Let a and h be the vectors parallel to the grain 
boundary segments ‘a’ and ‘b’ respectively. 
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From the figure 2.5, the vectors r, ,^2 ^ expressed as follows: 

rx=j 

^2 =-sina z +cosa j 
h~sm.9 cos (j)i +s\xv 6 sin ^ j + cos 6 k 

(2.9) 

since a is the line of intersection of planes A and P, while b is the line of intersection of 
planes B and P, a and b can be represented by the following cross-products: 

a = f^xn ( 2 . 10 ) 

y A A 

0 = ^2 X ^ 

( 2 . 11 ) 

from the equations 2.7 to 2.1 1, a and b can be expressed in the following form, 
a = cos6 i -sinOcoscj) k (2.12) 

-* A A 

b = cos a cos 6 i + sin a cos 9 j- sin 9 cos((zi -a)k 


(2.7) 


(2 8 ) 


(2.13) 
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15 , the plane dihedral angle, can be expressed as: 


cos = 


a • b 




a * 

b 


(2.14) 


The terms in above equation can be obtained from equations 2.12, 2.13 and 2.14: 


a-b = (cos(9 i -sin6cos<f) k)(cosacos0 i +sinQ;cos(9 j - sm.6 cos{(p - a) k) 


a -b = cos^ 0 cosa + s\n^ 6 cos (j) cos(^ - a) 


(2.15) 


\a\ 


-v/cos ^ ^ + sin ^ ^ cos ^ <(> = sjl - sin ^ 6* sin ^ (j) 


(2.16) 


= -Jcos^ acos^ ^ + sin^ crcos^ ^ + sin^ 9cos^{(p-a) = -yjl-sin^ 0 sin^(^-a) 


(2.17) 
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Thus from equations 2.8 to 2.1 1, the following expression can be obtained 


„ cos^ ^cosa + sin^ ^cos^cos(^-Qr) 
cos P = ■' ... Z X L. 

•y (1 - sin ^ (9 sin^ jzi)(l - sin^ ^ sin^ {(p - a)) 


(2.18) 


It IS seen from equation 2.12 that if ^ = 0 then P=a , as would be expected. However, 
for any given value of a , can vary from 0° to 180° depending upon the orientation of 
the plane of polish. 



z 



Figure 2.5; A random section through two grain boundary planes at an angle of a . 
The normal to the plane of section makes an angle 9 with the Z-axis and the 
projection of the normal on the XY-plane makes an angle (p with X-axis. The line 
of intersection of grain boundaries on plane of polish makes an angle . 
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On sectioning the polycrystalline structure (set X), the probability of an intersection of a 
triple edge (set Y) with the sectioning plane is given by (Sangal, 1998): 

P,(T 2 n YIT 2 tX) = — (2.19) 

71 


where / is the mean length of a triple edge. Consider a set consisting of pairs of 

intersecting planes making equal angles a distributed in the polycrystallme space. The 
intersection of Y with the sectioning plane, T 2 will give rise to a probability distribution 
of plane dihedral angles given by ! a)dfi . Let the numerical density of the set be 

'H\!{a)da , i.e. number of plane pairs making an angle of a per unit volume. The 
probability density, g( /? ) of plane dihedral angles is given by 


|P(T 2 n 7 1 Tj t Y) piP ! a)]<\y{a)da 

g(B)- 

Jp(T 2 nFITj t Y)Nv(a)^fa 

a 


( 2 . 20 ) 


Substituting P,(T 2 n Y | Tj t X) from equation 2.19, the probability density g{P) 
reduces to : 


g{l3)^ \p{fia)f{a)da (2.21) 

a 

Where f{a) = Ny{a)l ^Ny{a)da is the probability density function of TDA. Equation 

a 

2 . 21 , which is an integral equation, relates the PDA distribution to the TDA distribution. 
In order to solve equation 2.21, it is discretisized by considering discrete distributions of 
plane as well as dihedral angles. Both the distributions are divided into several classes. 
The integral equation 2.21 is converted into set of simultaneous equations: 
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Pxx Pa--- fx 

Px2 Pa--- fi 

g,=\ \ \ \ \ (2 22) 

■ Pxj f, 


til 

Where, g, is the frequency of the j class of plane dihedral angles; f is the frequency of 
the i^'’ class of true dihedral angles; p is the probability of generating plane dihedral angle 
in j'*’ class from a true dihedral angle in the i'*’ class. A Monte Carlo procedure has been 
developed for the estimation of coefficient py. For different class widths, a new set of 
values for these coefficients has to be calculated The frequencies gj can be determined by 
the measurement of PDA on the polycrystalline microstructures. Now the frequencies of 
TDA can be calculated by solving the simultaneous equations given m relation 2.22. 


2.1.3 Serial Sectioning and Disector Technique 

Test probes of fewer than three dimensions samples a particle in a microstructure with a 
chance proportional to their heights for a two dimensional probe, test lines intersect 
particles depending on their surfaces and test points samples particles proportional to 
particle volume. This method introduces a bias in the results and becomes an indirect way 
to count number of particles. Thus three-dimensional test probes were introduced to 
overcome the statistical errors caused due to the bias mentioned above. Measurement of 
average grain volume and certain topological parameters have been attempted by using 
three-dimensional serial section analysis by Rhines et al. (1976). 

Serial sections consist of consecutive slices cut through a piece of material at intervals 
that are so narrow such that every significant detail of internal structure are exposed as 
shown in figure 2.6. For topological measurements it is not necessary that the slices are 
parallel to each other, but for determining the total volume of material that has been 
sectioned it is essential that the sections are flat and parallel to each other. The total 
volume included in the serial sections, thus produced is calculated by making micrometer 
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measurements of the total thickness of the sample at the first and the last section. The 
number of features are counted in the sectioned material is then divided by the measured 
volume to obtain Ny (number of features per unit volume). Not only precision is to be 
maintained while making successive sections but also care is to be exercised while doing 
etching, repeated etching on a section can systematically alter the surface contour and 
may alter the quantitative observations. The counting of features on section is done by 
numbering each grain in polycrystals as shown in figure 2.7 Additionally it is 
recommended to do systematic bookkeeping to record the data in usable form. The 
amount of effort required for the entire procedure to be successful is enormous and fairly 
time consuming. 

By serial section method entire grain in a polycrystal can be reconstructed and true 
dihedral angles can also be evaluated. However this procedure was prone to many errors 
while tagging the series of sections accurately in an opaque sample This procedure 
turned out to be extremely time consuming. 



Figure 2.6 : computer animation of 
part serial sections stacked together. 



1 



IT 

\ V ;» 

\ &■ 



. 14 1 .5 / 


‘4 ,r^,. 




. jiv 

cs£i t » /■-L 


4 5 A 

Figure 2.7: Tracing of few representatives of serial 
sections. 
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Disector 


Sterio (1984) introduced the disector technique to overcome above limitations. In this 
technique two parallel 2-dimensional sections separated by a known thickness, h, are 
made on the sample and the counting is done as per requirement on these two sections 
The thickness is kept smaller than the smallest particle size in the microstructure The 
main advantage of disector is that nothing is assumed about the particles shape size or 
volume and actual counts are performed in the sample. Although the information 
obtained by this method is lower than serial sectioning but the overall efficiency 
increases multifolds, Gundersen (1986). The number of features appearing on one section 
and not on the other is counted which is divided by volume of the disector to give Nv 
(volume will be disector depth multiplied by the area of the section). This technique 
however demands accurate calculations of the disector depth, which is discussed in the 
following chapters. 


Dihedral angles 

Measurement of the geometry of triple junctions between grain boundaries in 
polycrystalline materials generates large sets of dihedral angles from which maps of the 
grain boundary energy may be extracted as shown in the figure 2.8, Adams et al. (1999). 
An extended form of orientation imaging microscopy (OIM) was used to measure both 
triple junction geometry via image analysis in the SEM and local grain orientation via 
electron back scatter diffraction. Serial sectioning with registry of both in-plane images 
and successive sections characterizes triple junction tangents from which true dihedral 
angles are calculated. The procedure gave an estimate of grain boundary energy but it 
turned out to be extremely tedious and inefficient due to serial sectioning but image 
processing introduced some ease in the evaluation. 
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Figure 2.8: Two planes of a polycrystalline microstructure 
indicating dihedral angles and grain boundary energy. (Murr, 
L.E). 


2.2 Digital Image Analysis 


The practical application of stereological methods depends upon the development of 
adequate technology. In 1916 Shand developed the first scanning stage for a light 
microscope in which the specimen was driven across the microscope field containing a 
hair line cross of two co-linear micrometers to measure lineal distances to measure 
volume fraction of a particular phase. Automatic methods of image analysis were 
introduced in quantitative analysis by recording a field under microscope with sensitive 
photo detectors. These detectors record image on the basis of contrast differences. A 
more powerful variant called the flying spot microscope was developed by Roberts in 
1951. Here a fine spot of light produced by a cathode ray oscilloscope was focused on 
the specimen and the image was rapidly scanned. A photomultiplier picked up the signal 
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emitted from the specimen and transferred to second cathode ray tube which was 
synchronized with the first tube. 

In 1967, TV cameras were introduced to cut up the microscope image into a dense lattice 
of fine lines or points. Computers and image processing algorithms hence provided 
powerful tools to examine microstructures. It pnmarily acts as tool to automate the 
procedure of determining the evaluation of tedious stereological parameters like dihedral 
angles. Digital image processing have been used by many investigators in the field of 
quantitative characterization of multi phase materials like Silva et al. (1998) and Latala et 
al. (2001). They have concluded that the results obtained by computer-aided analysis 
provide highly consistent results. Though computer aided processing takes somewhat 
longer time to set up various parameters but then due to fully automated further analysis 
the total analysis period was about 100 times shorter than manual computations. Gokhale 
et al (1999), developed computer based techniques to determine grain size distribution as 
well as grain volume determination using serial sectioning and large area disector probes. 


2.2.1 Image Digitization 

An image is visual representation of an object, and image processing manipulates 
information within an image to make it more useful. In order to process an image with a 
computer the image must be converted into numeric form, this process is termed as image 
digitization. The digitization process divides an image into horizontal grid of very small 
regions called picture elements or pixels. In computer thus the whole image is 
represented by this digital grid. This grid is more generally called as ‘bitmap’. Each pixel 
in the bitmap is identified by its position in the grid by it’s row (x) and column number 
(y). The convention used (in this study) is that the (0,0) pixel is located at the upper left 
position of the bitmap as shown in the figure 2.9. when the source sample is scanned 
under a CCD camera it is examined in a grid fashion. That is each pixel in the image is 
individually sampled and the brightness is measured and qualified. This measurement 
results in a value for a pixel, for a gray scale 8 image, it will range from 0 - 255. Gray 
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scale-8 denotes that 8 bits are used to store each pixel value, which is also termed as the 
pixel depth. Depending upon the number of bits used for storing pixel value, the image 
can be classified as Gray Scale 8, Gray Scale 12 or Gray Scale 16. a totally black and 
white image is represented by 1 bit only, either the pixel value is ‘0’ completely black or 
the pixel value is ‘ 1’ completely white. In an 8 bit Gray scale image, a pixel with value of 
zero is completely black and a pixel with a value of 255 is completely white. A value of 
127 represents a gray color halfway between black and white (medium-gray), and a pixel 
value of 64 has a color halfway between medium gray and black. 


Pixel (0,0) 



Figure 2.9; A digitized image indicating array of pixels, 
the pixels in the above microstructure are shown much 
larger than the actual size. 
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2.2.2 Image Enhancement 


The acquired image requires some form of enhancement to improve its appearance and to 
extract useful data from it. Enhancement technique range from simple operations such as 
brightness and contrast adjustment to more complicated such as spatial and 
morphological filtering operations. These operations are able to refine the visual 
information as per the need. The image enhancement techniques are broadly classified 
under following categories: 

a) Modify the intensity values of the pixels, or Point Operations 

b) Application of spatial filters 

c) Manipulate the image frequencies 

Above enhancement process can be applied individually or in combo to achieve the 
desired results. 


2. 2. 2.1 Point Operations 


The image can be enhanced to a limited extent by modifying the intensity values of the 
pixels. If the image is very dark overall, all values of the pixels can be increased by 
stipulated amount to make it better. Brightness, contrast and thresholding are primary 
operations, which work by altering the intensity values of the pixels. These 
transformations are termed as point operations as the result of any transformation of any 
pixel depends only on the initial gray value of this pixel and is independent of it’s 
neighbors. When brightness is increased the value of each pixel of the image is increased 
closer to value ‘255’ in a Gray Scale 8 image. When brightness is decreased the pixel 
intensities are reduced, moving them closer to 0, or black. 
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Contrast denotes the degree of difference between the brightest and darkest components 
of an image. An image has poor contrast if it contains only harsh black and white 
transitions or contains pixel values within a narrow range. An image has good contrast if 
it is composed of a wide range of brightness values from black to white. Increasing the 
contrast spreads the pixel values across a wider range, while decreasing contrast 
compresses the values into a narrower range. 


Thresholding allows an image to be converted into just two colors: black and white. This 
is done by specifying a range of intensities to be emphasized, they are set to white and 
converting all others to black thresholding is used for segmenting an image to extract 
important features for further measurements or to reduce an image to two intensity levels 
m preparation for filtering operations. 


Binarization 

Currently gray images are most frequently used aid for recording image data and 
processing it. Binarization is an important step in image analysis (Wojnar, 1998) as the 
binary images occupies approximately 8 times smaller memory than gray images and 
almost 24 times shorter than a colored image also only in binary images can one detect 
separate features such as grains (every connected set of pixels is recognized as a single 
particle). These separated objects can now be subjected to measurements such as area, 
perimeter, diameter etc. another advantage of binary image is that some transformations 
mainly morphological operations can de done on such images only. There are procedures 
to separate connected particles, to reduce grain boundaries to one pixel thin lines etc 
falling under this category. The process of transformation of gray scale images to binary 
images is called binarization or thresholding as shown in figure 2.10. 
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2.2.2.2 Spatial Filtering 

Filtering operations reduce or increase the rate of change that occurs in the intensity 
transitions within an image. Areas, in which there are sudden changes in intensity, appear 
as distinct edges in the image. Area where there are gradual changes shows as soft edges 
Filtering acts to detect and modify the rate of change at these edges. It can increase the 
intensity difference s in a soft image to make it sharper or reduce the intensity differences 
in a hard edge to smooth and soften it. 

Filtering operations produce their effect by modifying a pixel’s value based upon the 
values of pixel’s surrounding it. For example, blurring is accomplished by averaging all 
of pixel values and replacing the center pixel with the averaged value. This reduces the 
variations among neighboring pixels. Filtering techniques axe divided into two categories: 
convolution and nonconvolution filters. Both techniques accomplish their results by 
examining and processing the image in small regions called pixel neighborhoods A 
neighborhood is a square region of image pixels, typically 3X3, 5X5 or 7X7 in size. 
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Convolution filters 


Enhancement or convolution filters process image neighborhoods by multiplying the 
values within a neighborhood by a matrix of filtering coefficients. This matrix is called a 
kernel. It is the same size as the neighborhood that it is being applied to. The results of 
this multiplication are summed and divided by the sum of the filter kernel. The result 
replaces the center pixel (sometimes also called origin), in the image neighborhood. The 
convolution process always uses a neighborhoods original pixel values as input. For 
example the kernel focus is moved to say, pixel 6 (on the image), the filtering process 
will use pixel 5’s original value, not the one it was just assigned by convolution. The 
various filters under this category are low-pass, hi-pass, sharpen median etc. 


Non-Convolution/ Morphological filters 

Morphological filters are highly abstract theory of image transformations, possessing it’s 
own rules and notations (Wojnar, 1998). Morphological processing refers to certain 
operations where an object is hit with a structuring element (kernel) and thereby reduced 
to a more revealing shape. The stmcturing element are defined using following notation: 

1- represents pixels belonging to the set of points under analysis for example black 
points in figure 2.1 1(b). 

0- represents pixels belonging to the matrix for example white points in figure 
2.11(b) 

X- represents pixels not taken into account, that is, this point can have any value and 
has no effect in the transformation. 

Non-convolution filters also work with pixel neighborhoods but do not multiply the 
neighborhood values by a kernel like in convolution filter. This filter works with the data 
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in the neighborhood itself and uses either a statistical method or a mathematic formula to 
modify the pixel upon which it is focused. Erosion , dilation, open and close are few 
examples of morphological filters, which is described in detail in section (experimental 
procedure). 
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CHAPTER 3 

MATERIALS AND MICROSCOPY 


3.1 MATERIAL 


The experiments were performed on commercial polycrystalline 316L austenitic 
stainless steel. The chemical composition of the steel (determined using JEOL 
electron probe micro analyzer) is given in Table l.The material, in the form of 
12.5 mm thick plates, was cold rolled to 92% reduction. 


Table 3.1: Nominal chemical composition of 316L austenitic stainless steel. 


Elements 

Si 

Ni 

Ti 

W 

C 

Cr 

Mo 

Cu 

A1 

Fe 

Composition 

% 

0.528 

13.73 

0.029 

0.076 

0.003 

18.05 

2 025 

0.019 

0.036 

balance 
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3.2 HEAT TREATMENT 

The cold rolled material was given five different annealing treatments as given in 
the Table 2. The heat treatment was carried out in a Kanthol furnace. Small pieces 
of rolled samples were annealed at five different temperatures. The Table also 
gives the variation in grain size (measured as mean intercept length). The method 
of determining grain size is discussed in section 2.1.1. 


Table 3 2: Statistical parameters of plane dihedral angles and annealing 

temperatures 


ANNEALING 

ANNEALING 

GRAIN 

DISECTOR 

TEMP 

TIME 

SIZE* 

THICKNESS 

(•c) 

(Minutes) 

{fjm) 

h^jLon) 

900 

60 

24.916 

3.80 

900 

180 

26.69 

4.50 

950 

60 

35.93 

3.18 

1000 

60 

53.39 

6.16 

1050 

60 

62 29 

6.25 


mean intercept length has been taken 
as a measure of grain size 
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3.3 Microscopy 


3.3.1 Specimen Preparation and Examination 


For accurate automated measurements of microstructural parameters sample 
preparation becomes critical. Mechanically polished samples were 
electrochemically etched: 60% HNO3 was used as electrolyte; electrode used was 
316L stainless steel; current density was kept at 0.75 mA/mm^. This etching 
procedure ensures that the twin boundaries are not revealed. Another advantage is 
that it preferentially etches the grain boundaries resulting in a sharp black and 
white contrast. Such a contrast is ideal for automatic image analysis. 


3.3.2 Image Analysis and Disector Thickness 

24 bit RGB digital images were acquired through a CCD camera mounted on 
Lieca Optical Microscope. These images were subsequently converted to 8-bit 
grayscale images. The size of the acquired image was 640X480 picture elements. 
A micro-scale was used to perform calibration of the image. At a magnification of 
lOOX, dimension of one pixel was calculated to be 0.6041 ^ . The software used 
was Image-Pro Plus ver 4.1 for capturing as well as processing the images. 


After polishing and etching the first section of the disector, several micro hardness 
indents were put as reference markers using a load of 160 gms by a Leica micro 
hardness indenter. The indentor used was a standard diamond pyramid having an 
angle of 136° between its two opposite faces. Digital image processing was 
employed in order to accurately measure microstructural parameters and the 
positions and sizes of the hardness indents. This was followed by mechanical re- 
polishing and re-etching to produce the second section (parallel to the first 
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section) of the disector. The hardness indentations were used to align the images 
of the two sections as shown in figure 3.1. The decrease in the size of the indents 
was used to establish the separation distance, h, of the two sections Si and S 2 as 
shown in figure 1 . The distance of separation, h, which is termed as the disector 
thickness, is given by : 

h = (D, - D 2 )/ 7.0006 (3. 3.2.1) 

where, Di and D 2 are the average diagonal lengths of indents in sections Si and S 2 
respectively 



Figure 3.1: The disector consisting of the two parallel 
sections Si and S 2 separated by a distance h. Indents are 
produced by a standard diamond pyramid indent of a 
microhardness tester Di and D 2 are the diagonal lengths of 
the indents. 
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CHAPTER 4 

IMAGE ANALYSIS ALGORITHM AND THE DISECTOR 
TECHNIQUE 


4.1 Image processing algorithm 

4.1.1 Binarization/ Thresholding 

A typical grayscale image of the microstructure of 3 1 6L austenitic stainless steel 
and the corresponding bimodal gray level distribution are shown in figures 4.1(a) 
and 4.2(b) respectively. For the binarization of such images as mentioned in 
section 2.2, automatic thresholding was applied, which involved fixing the 
threshold limit at a local minimum between two local maxima of gray level 
distribution. All pixel values below the threshold limit were set to ‘0’ value and 
those pixels whose value was above this limit were fixed to ‘255’ value. This type 
of binarisation gave consistent results independent of image contrast and 
brightness. Figure 4.1(c) shows the resulting binary image in which all the sensor 
noise and background illumination variation have been abstracted. 
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(b) 



(c) 


Figure 4. 1 : (a) A typical grayscale image of a polyciystalline microstructure; (b) 
Bimodal gray level histogram of the image in figure (a), arrow mark indicates 
thresholding limit; (c) transformed binary image (after applying thresholding 
algorithm) showing polishing/etching artifacts 
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4.1.2 Erosion filtering 


The erosion filter is a morphological filter that changes the shape of an object in 
an image by eroding or reducing the boundaries of the specified objects. Suppose 
the object X and the structuring element B are represented as sets in Euclidean 
space. Let Bx represent translation of B so that it’s origin is located at x. then the 
erosion of X by B is defined as the set of all points x such that Bx is included in X: 

Xi-)B A {x-.B^dX} (3.1) 


Erosion is a shrinking filter therefore eroded objects are smaller than the initial 
one. Some narrow peninsulas or small particles will disappear after erosion. 
Figure 4.2 shows example of erosion filter. 


This filter removes noise/artifacts created by etch pits, inhomogeneous 
etching/polishing etc, in the image, erosion filter was applied which removed all 
the pixels touching the matrix. This operation resulted in the removal of some 
narrow peninsulas and small artifacts as shown in figure 4.3. A typical erosion 
kernel is shown in figure 4.3(a), this kernel or the structuring element fits with the 
configurations on the image pixels and does transformation only to the pixels 
which matches exactly with this structuring element. It is often used to reduce , or 
eliminate small bright objects as this filters erodes the boundaries of bright objects 
and enlarges the boundary of the dark ones. The image shown in figure 4.3b has 
been inverted to make convenient comparison with the original image as shown in 
figure 4. 1 (a). 


Object X 


Structure 

element 


Eroded object 
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Figure 4,2: illustration of erosion of object X by the given structural 
element, notice the blank circles indicating erosion of the boundary 
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4. 1.3 Thinning Morphologicai Filter 


The eroded image is transformed by a thinning morphological operation to reduce 
the grain boundaries to only one pixel thick line. The thinning filter was applied 
Iteratively till further operations of the filter had no effect on the image. The 
filter’s basic operation is as shown in figure 4.5. The resulting thinned image is 
shown in figure 4.6. However, this image had numerous open-ended branches, 
which appear because of non-homogeneous etching Such branches can lead to an 
erroneous detection of triple points. 





Structuring 

element 


0 ^ ^ngm 

• •• 




Thinned 


• OCClOCClOC*COO» 

Cl*'3COw'CCiOOOO«0 

* 0 0 0 c 0 c 0 0 0 0 i 


Figure 4.5: Illustration of thinning morphological transformation of the 
object, notice the reduction of the object to one unit thick line. 



Figure 4.6: Image of microstructure resulted after thinning 
morphological filter was applied on image shown in figure 
3.3, notice that the grain boundaries now are only one pixel 
thick. 
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4. 1.4 Pruning Filter 

Therefore to smoothen the thinned image and to remove all parasitic branches a 
pruning filter was applied. Figures 4.7 and 4.8 show the basic pruning operation 
on the object and the final pruned image respectively. 



Figure 4.7: Illustration of pruning transformation of the object, notice that the 
parasitic branch has disappeared. 



Figure 4.8: Image of microstructure resulted after pruning 
filter was applied on image shown in figure 4.6, notice that 
the grain boundaries now are only one pixel thick. 
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4,1.5 Bitmaps 


The final image was converted into a bitmap. Figure 4.9 shows a bitmap of the 
rectangular region marked in figure 4.8. As can be seen fi-om this figure, the grain 
boundary pixels are represented by O’s and the matrix pixels are represented by 
1 ’s (these pixels uniquely had a value of 255 in the binary image). The bitmap was 
used as input to the algorithm developed to detect triple junctions and determine 
their geometry, as discussed below. 


111111111101111111111111111111111111111101111111111111111111111111111111 

111111111101111111111111111111111111111101111111111111111111111111111111 


111111111 

111111111 


roixmiiTiJiiiiiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiii 111 111 11 1111111 

110111111111111111111111111111101111111111111111111111111111111 


111111111111000001111111111111111111111101111111111111111111111111111111 
111111111111011100000111111111111111111101111111111111111111111111111111 
iiiiiiiiijiiojLniiiiiJioooiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiiiiiiiiiiiiiiiiii 
iiiiiiiinoinrrriniiiioooiiiiiiiiiiiiioiiiiiiiiiiiiiiii 11111111 1111111 
111111111011111111111111111000111111100011111111111111111111111111111111 
111111111011111111111111111111000001011101111111111111111111111111111111 
111111111011111111111111111111111110111101111111111111111111111111111111 
111111111011111111111111111111111110111101111111111111111111111111111111 
111111110111111111111111111111111111011101111111111111111111111111111111 
111111110111111111111111111111111111011101111111111111111111111111111111 
111111110111111111111111111111111111100010111111111111111111111111111111 
111111101111111111111111111111111111101111011111111111111111111111111111 
111111101111111111111111111111111111101111011111111111111111111111111111 
111111101111111111111111111111111111011111101111111111111111111111111111 
111111101111111111111111111111111111011111110111111111111111111111111111 
111110010111111111111111111111111111011111110111111111111111111111111111 
110001111001111111111111111111111110111111111011111111111111111111111111 
001111111110011111111111111111111101111111111011111111111111111111111111 
111111111111100111111111111111111011111111111101111111111111111111111111 
111111111111111001111111111111110111111111111101111111111111111111111111 
111111111111111110011111111111101111111111111110111111111111111111111111 
111111111111111111101111111111011111111111111110111111111111111111111111 
111111111111111111110011111110111111111111111111011111111111111111111111 


Figure 4.9: Bitmap of the marked rectangular region in 
figure 4.8. ‘0’ refers to grain boundary pixels and ‘ 1 ’ 
refers to matrix pixels. 
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4.2 Image Analysis Algorithm: Evaluation of Plane 
Dihedral Angles 


4.2.1 Triple Point Detection 


Refer to the bitmap of the rectangular section of the image (figure 4.8) of 
microstructure as seen in fig 4.9. Each grain boundaiy pixel (pixel value of ‘0’) 
was considered at the center of a 3x3 matrix. If a grain boundary pixel had exactly 
three other boundary pixels in its neighbourhood then the pixel was marked as a 
possible triple point and subjected to further evaluation by comparing it with a 
pre-generated list of valid triple point configurations. A partial list of valid 
configurations in the form of 3x3 matrices are shown in figure 4.10. The 
remaining configurations can be generated by rotating the individual matrices 

7t 'hK 

through -^,71 and — thereby generating a total of 12 valid configurations. 

However there are many other configurations possible in which the central ‘0’ 
pixel has three neighbouring pixels, but are not included in the above list. Some 
examples are listed in figure 4.11. The reason for this becomes evident from 
figure 4.11, which shows a magnified bitmap view of the rectangular region 
marked in figure 4.9. The (row, column) coordinates of each pixel in the bitmap 
can be referred relative to the left-top pixel, whose coordinate is assigned (1,1), as 
shown in figure 4.11. On examining this figure, it is clear that grain boundary 
pixels at coordinates (3,6), (3,7), (4,7) and (4,8) are not triple points even though 
each of these pixels have exactly three neighboring grain boundary pixels, such 
triple points have been termed as pseudo triple points. Similarly the grain 
boundary pixels at coordinates (3,4) and (4,3) are also pseudo triple points. 
However the triple point at coordinate (3,3) is an actual triple point, which 
matches with the configuration shown in figure 4.10a rotated through n: . Thus, 
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the pseudo triple points can be classified according to the following general rule: 
If any two of the three neighboring pixels share a boundary among themselves, 
which means that if any two neighbouring pixels have same ‘x ’ coordinate and the 
‘y ’ coordinate differs by one and vice-versa then the triple point is classified as a 
pseudo triple point This rule was found to be appropriate by comparing the 
results of the automated detection of triple points with manual detection. 



Figure 4.10. A partial configurations of valid triple points. Remaining 
configurations can be generated rotating the matrices through tv 12, 71 : and 'hit! 2 
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Figure 4.11: Illustration of pseudo triple points (refers to 
rectangular region marked in figure 4.9) 
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4.2.2 Evaluation of Plane Dihedral Angles (PDA) 


Starting from each of the detected triple points, as described in section 4.2.1, the 
coordinates of the pixels of the grain boundary segments were traced till the 
segment terminated into another triple point or till no further touching pixels could 
be found. The algorithm developed for the same is at appendix A.l. At this stage, 
if a particular segment was found to be very short (less than four pixel length) 
then the associated triple point was considered as a pseudo-triple point and 
therefore eliminated from the analysis. Manual examination of the micro structure 
confirmed the validity of this procedure. The output of the program is at appendix 
B.l. 

A least squares procedure was employed to fit polynomials on the pixels of the 
grain boundary segments at each triple point. The form of the polynomial chosen 
for the purpose is given as: 


y = yo + a^ix - x^) + a^ix - Xa)^ +...-t-a„(x-Xo)'” (4.2.2.1) 

Where ‘m’ is the polynomial order, ai, a 2 am are the coefficients. This form of the 
polynomial was modified to ensure that the fitted curve is constrained to pass 
through the triple point. 

A least squares procedure was adopted to determine the values of the polynomial 
coefficients. This involved minimization of the sum of squares (R^) of the 
following equation: 


(=1 


(4.2.2.2) 


where, y, is the measured value and y is given by equation (1) for the ith pixel. 
Derivative of equation 4.2.2.2 with respect to each of the unknown coefficients of 
the polynomial is performed and equated to zero to develop a set of equations 
which was solved by gauss-elimination technique. 
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Each grain boundary segment was fitted with polynomial of increasing degree of 
order, starting with the order of 1. In order to select the best polynomial order, a 
ratio F was evaluated to compare the goodness-of-fit (quantified by sum of 
squares, R^) of two polynomials of orders m and m+1 


F = 


{Rl-Rl)/ 

(ri -Yi)/ 
/Yx 


(4.2.2.3a) 


and 


y, = « - (m + 1) and =n-(m + 2) 


(4.2.2.3b) 


where, y = degree of freedom, n = number of data points and the subscripts 1 and 2 
refer to the polynomials of orders m and m+1 respectively. 

If the fit of the higher order polynomial (order = m+1) has a larger sum of squares, 
i.e., R 2 ^>Ri^, then clearly the fit becomes worse on increasing the order and 
therefore, the lower order polynomial (order = m) is chosen as the best order. In 
terms of the F-ratio (equation 4.2.2.3a), the value of F will be negative. However, 
in general with increasing polynomial order the sum of squares value is expected 
to reduce because the higher order equation has more adjustable parameters. Thus, 
a further analysis of the F-ratio was done. 

In going from lower order to the higher order polynomial, if the relative decrease 
in the sum-of-squares was less than the relative decrease in the degree of freedom 
(i.e., F<=1), then the polynomial with the lower order was selected as the best 
polynomial order. On the other hand, if the value of F was greater than 1, i.e., the 
relative decrease in the sum of squares was larger than the relative decrease in the 
degree of freedom, then the higher polynomial order was accepted and the above 
analysis was repeated by further incrementing the polynomial order. 
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The correlation coefficient (Chapra 1989), r, for each polynomial fit was 
determined by using the following equation; 


r 


2 


S^-R^ 


( 4 . 2 . 2 . 4 ) 


where, =2(T, - y? ■ 

i=l 


For a perfect fit = 0 and = 1, signifying that the poljmomial fit truly fits all 
the coordinates of each grain boundary. While r^ close to zero signifies a poor fit. 
The correlation coefficient was used to decide whether to accept or reject a 
particular fit, as discussed in the section (chapter 5 results and discussions). 


The plane dihedral angles were determined as the angles between the tangents to 
the grain boundary segments at the triple points. Figure 4.12 illustrates a typical 
example of three tangent lines (calculated by differentiating the fitted curves) at a 
triple point. Plane dihedral angle, , is given by 


cos P = 


(OAf +(OBy -(ABp 
2{OA)(OB) 


( 4 . 2 . 2 . 5 ) 


Similarly, the other two plane dihedral angles can also be determined. 
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Figure 4.12: Tangents (thick lines) drawn to the fitted grain 
boundary segments (thin lines) at a triple point ‘O’. ‘A’ and ‘B’ are 
arbitrary points on the tangents. 
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4.3 3-Dimensional Stereological Probe: The Disector 


4.3.1 Transformation of Plane Dihedral Angles (PDA) to True 
Dihedral Angles (TDA) 


Figure 4.13 shows the geometry of a disector cutting the grain boundary planes 
meeting at a triple edge. The two sections, Si and S 2 (with a separating distance of h) 
are parallel to the x-y plane. The vectors and are the unit normals to the grain 
boundary planes 1, 2 and 3 respectively. The orientation of the triple edge is defined 
by the vector.^ . The intersection of the disector with the grain boundary planes result 
in linear grain boundary segments meeting at triple points ‘Q’ and ‘P’ on the parallel 
sections Si and S 2 respectively. The orientation of the grain boundary segments PA, 
PB and PC are defined by the unit vectors and respectively. The plane 

dihedral angles, and , are the angles between the vectors r^ and and 
fj, and Fjand fy respectively. The coordinates of the triple points ‘Q’ and ‘P’ are 
denoted as (txi.tyi) and (tx 2 ,ty 2 ) respectively. 
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Section S i 



Figure 4.13: Random disector section of thickness h meeting at a triple edge PQ. 
Coordinates of P (tx 2 , ty 2 ) and coordinates of Q: (txi, tyi) 
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The true dihedral angles and are given by the following general equations 

(refer to figure 4. 1 3): 


. . £ X f. 

cosa, =«,•«,= — i- 

£xfi 


£xf, 

£x fj 


cos (22 = ■ ”3 


£ X f, £ X n 


£xrl 


£xn 


(4 3.1.2a) 


(4.3.1.2b) 


a-^-lK- (a, + (Xj) 


(4 3.1 2c) 


The following transformation of coordinates can be applied without any loss of 
generality (see figure 4.13). The origin is translated from point O to the triple point P 
followed by a rotation about the z-axis (angle between the grain boundary segment 
PA and the x-axis), through ^ which is given by 


coscj) = 


(^1-^x2) 


[(^1 -hif -^iy I -tyif} 


or 


sm 0 = ■ 


iy^ -tyl) 


[(Xj t^2) + (Ti ^yl^ ] 


(4.3.1.3a) 


(4.3.1.3b) 


where, Xi and yi are the coordinates of an arbitrary point on the segment PA. The 
above transformation results in the unit vector r^ coinciding with the transformed 
x^-axis. In this transformed system of axes (x^, y, ?! ) the coordinates of the triple point 
P is (0,0,0) and the coordinates of triple point Q is (tx, ty, h): 



h = ihx - 1.2 ) COS ) sin (!) 


(4.3.1.4a) 


ty =(^i -^y2)cos<zi-(/,, -t^^)sm<l) 


(4.3.1.4b) 


Thus, from figure 4.13, the vector E and the unit vectors f, and are given by the 
following expressions: 


E = it, + jty + kh 

(4.3.1.5a) 

^ • (4.3.1.5b) 

Fj = i cos + y sin /?, (4.3.1.5c) 


The true dihedral angle can be obtained by substituting equations (4.3. 1.5) in 
equation (4.3.1.2a): 


cos a, = 


jh - kt 


z /z sin + jh cos + k{t^ sin 1 cos y^, ) 




[h +(r^sin^i-t cos>9,)“] 


2il/2 


On simplification : 


cos a, = 


if cos y5j - ty sin y9, - ty COS y9, ) 
{h^ +tiy'^ [If + (t^ sin y9, - 1 cos 


2 -i \/2 


(4.3. 1.6) 
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Similarly equations fora 2 can be derived by performing a transformation of 
coordinates in which the rotation about the z-axis is done through (p + (3^^ (instead of 
(j ) , as above) 


The expression for 0^2 (which has the same form as equation 4.3. 1.6) is given below: 

cosa cos /; (/; sin - 1'^ cos fS-, ) 

{h^ +(t[sm/3^-t'yCosP^ff^ (4.3. 1.7) 


where t[ andf^ (similar to tx and ty in equations 4.3.1 .4) are: 


t'x =(^i -^. 2 )cos((^ + /?i) + (r^, -t^ 2 )sin((z) + ;0i) (4.3.1.8a) 

t'y = (^1 - ty 2 ) cos{(j) + A ) - - 1 . 2 ) sin(<zi + A ) rA■x^ RK^ 


Above derived disector technique has been used to estimate true dihedral angles from 
digital images, as discussed below. 


4.3.2 Computation Methodology of True Dihedral Angle 


Figure 4.14a shows a typical unprocessed grayscale image of a polycrystalline 
microstructure of 316L austenitic stainless steel annealed at 900° C for 60 minutes. 
This image is the first disector section Si as already explained in figure 3.1. The 
figure also shows the micro hardness indents which were used as reference markers. 
The grayscale images were binarized, etching/polishing artifacts were removed, and 

frnRfiPzT slratrTTrT 
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the grain boundaries were thinned and pruned by an image-processing scheme 
developed in section 4.1. The result of the image-processing scheme on the image in 
figure 4.14a is shown in figure 4.14b. The sample was re-polished and re-etched to 
generate the disector section S 2 , as shown in figures 4.14c and 4.14d (processed 
image). The size of the indents in sections Si and S 2 (see the unprocessed images of 
the two sections in figure 4.14a and 4.14c respectively) were substituted in equation 
3.1 to determine the disector thickness, h. For accurate measurement of the disector 
thickness the indent size in section S 2 was measured after polishing but prior to 
etching. This avoided any errors resulting due to the distortion of the indent due to 
etching. The various disector thicknesses obtained for samples subjected to different 
annealing treatments are shown in table 3.2. 

In order to match the center of the corresponding indents in the two images of 
sections Si and S 2 with each other precisely, an annotation mark was made at the 
center of the indents. The images of the two sections were merged such that the 
centers of the indents were matched. Care was taken to mount the sample in such a 
way that almost no rotation was needed to match the center of the indents. The result 
of the match of the first and the second image was found to be within 2 to 3 pixel 
accuracy, as shown in the composite image in figure 4.14e. The light image is from 
section Si and the dark image is from section S 2 in figure 4.14e. Triple point pairs 
similar to the triple point pair P and Q shown in figure 4.13 is shown in figure 4.14f 
which is enlarged image of the rectangular section marked in figure 4.14e. For 
example the pair of triple point marked as ‘r(belongs to Si) and ‘2’ (belongs to S 2 ) 
are from a single triple edge. The varying relative shift (in the 2-dimensional plane) of 
triple points in section S 2 , with reference to these on section Si is done to the random 
orientation of the triple edges. It may also be noted that there are many triple points in 
section Si (see rectangular region in figure 4.14f), which do not have any 
corresponding point in section S 2 . This is because of improper etching in section S 2 it 
could also arise if the disector thickness is larger than the size of the grains. 

The measurement of the coordinates of triple points orientation of the grain boundary 
segments, and the plane dihedral angles in the binarized images of plane sections Si 
and S 2 were performed by an automated image analysis algorithm described in section 
4.2. Corresponding triple point pairs in the two sections were identified by an 
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algorithm, which used the displacement of triple points in section S2 relative to these 
in section SI, and the closest neighbors were paired (program in appendix A.5). This 
simple algorithm was successful in microstructural images containing large grains 
(such as the structure obtained in samples annealed at temperatures of 1000° C and 
above). For microstructures with smaller grains the procedure was only partially 
successful and manual intervention was required. This problem could be eliminated 
by additionally comparing the orientation of the corresponding grain boundary 
segments. However, no further attempts were made in this direction. Finally, 
equations 4. 3. 1.6, 4.3. 1.7 and 4.3.1.2a were used to estimate the true dihedral angles 
and their distributions. 
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Figure 4.14. Microstxuctures of 316L austenitic stainless steel with microhardness indents as 
reference markers: (a) unprocessed image of disector section Si. (b) Processed image 
of disector section Si; (c) Unprocessed image of disector section S 2 ; notice the reduction in the 
sizes of the indents;(d) Processed image of section S 2 . (e) Merged sections SI and S2 using cenU 
of indents as reference; (f) image of marked section in figure 2e indicating a random triple edge. 
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CHAPTER 5 

RESULTS AND DISCUSSIONS 


As mentioned in chapter 3, the commercially available 316L austenitic stainless steel was 
obtained and were cold rolled to reduce thickness by 92% in more than 25 successive 
passes. To achieve range of grain sizes the cold rolled samples were annealed at different 
temperatures and times as indicated in table 5.1. The microstructure of the samples was 
characterized after capturing digital images and processing them as shown in figure 5.1. 
In general the digital images were converted into bitmaps and was used as input to the 
computer algorithms developed to evaluate dihedral angles. The programs are listed in 
appendix A and the outputs of these programs are illustrated in appendix B. 


In order to ensure measurement accuracy of plane dihedral angles by the automated 
technique, the correlation coefficients of the best-fit polynomials were checked. Figure 
5.2, shows typical example of the polynomial fit produced along with the correlation 
coefficients at a triple point. A statistical sample of manual measurements was compared 
with automated measurements along with the correlation coefficients (figure 5.2). By 
such comparisons it was observed that in the case of each large discrepancy between 
manual and automated measurement the correlation coefficient for one or more of the 
best-fit polynomial was less than 0.40. Thus for all subsequent automated measurements 
of plane dihedral angles, if the correlation coefficient was less than 0.40 then the 
associated triple point was rejected from the analysis. In the analysis of the 
microstructures of the samples of 316L austenitic stainless steel such a procedure resulted 
in approximately 5% rejection of the triple points. Thus by using the above technique, it 
was possible to obtain a statistically significant sample of the plane dihedral angles. Table 
1 gives the total number of triple points identified and the valid triple points finally 
evaluated per image. 



56 


Considering all the measured plane dihedral angles, histograms of the relative frequencies 
of plane dihedral angles were generated. Figure 5.3 shows a histogram of the relative 
frequencies of plane dihedral angles for the sample annealed at various temperatures 
Figure 5.4 shows the variation of standard deviations of plane dihedral angles as a 
function of annealing treatment. From this figure it is evident that with increasing 
annealing temperature, the spread of the distribution reduces about the mean value of 
120°. This is expected because with increasing annealing temperature grain boundaries 
tend to equilibrate [4]. 


Table 5.1: Statistical parameters of plane dihedral angles 


NN 

BMP. 

■c) 

ANN 

TIME 

(Minutes) 

GRAIN 

SIZE* 

i/um) 

STATISTICAL PARAMETERS 

NO OF 

TRIPLE 

POINTS 

EVALUATED 

STANDARD 

DEVIATION 

(PDA) 

DISECTOR 

THICKNESS 

ijum) 

STANDARD 

DEV 

(TDA) 

STD. 

DEV 

(GBE) 

MEAN 

(GBE) 

900 

60 

24.916 

199 

26.0393 

3.80 

26.95 

0.186 

0.57 

900 

180 

26.69 

160 

26.508 

4.50 

24.64 

0.187 

0.570i 

950 

60 

35.93 

110 

24.89 

3.18 

23.38 

0.151 

0.62 

1000 

60 

53.39 

120 

24.72 

6.16 

18.79 

0.14 

0.70 

1050 

60 

62.29 

52 

23.2 

6.25 

18.153 

0.139 

0.744 


* mean intercept length has been taken as a measure of grain size 
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As discussed in section 4.3 the two disector sections Si and were processed as 
mentioned in section 4.2 to obtain plane dihedral angles. Also, Corresponding triple point 
pairs in the two sections were identified by an algorithm (appendix A.5), which used the 
displacement ol triple points in section S2 relative to those in section SI, and the closest 
neighbors were paired to form a grain boundary edge as shown in figure 4.13. Equations 
4.3. 1.6, 4.3. 1.7 and 4.3. 1.2 were used to estimate the true dihedral angles and their 
distiibutions. Histograms of the relative frequencies of the true dihedral angles were 
generated, which is shown in figure 5.5. Figure 5.6 shows the variation in standard 
deviations as a function of annealing treatment. From the figures it was evident that as the 
annealing temperatures were increased the standard deviations decreased and the peaks 
shifted towards 1 20 class. 1 his explains that at higher temperatures the grain boundaries 
tend to equilibrate as also mentioned in the case of plane dihedral angles. Thus in the fine 
grain microstructures the standard deviations of TDA is higher and the distribution peak 
is away from 120'’ class in comparison to coarser grains microstructure. This indicates 
that in the fine grain structures there is a state of non-equilibrium at the grain boundaries. 

On comparing the standard deviations of the plane dihedral angles with corresponding 
standard deviations of true dihedral angles it was found that in general (particularly 
samples annealed at lOOO” and 1050” C) the standard deviation of TDA was lower than 
corresponding standard deviation of PDA. This can be understood in terms of the factors 
contributing to the spread of the distribution. The standard deviation (or variance) of the 
plane dihedral angle distribution is affected by the variation in the true dihedral angles 
and the variation due to sampling of randomly oriented triple edges by the sectioning 
plane (see Figure 2.4b). 'fhe latter contribution is eliminated during the process of 
transformation from plane to true dihedral angles leading to a relatively smaller spread in 
the distribution, 'fhis fact would become significant in studies attempting to relate the 
factors, such as, deformation temperature, strain, etc. to the distribution of dihedral 
angles. 

The samples annealed at 900” C, the standard deviation of TDA is comparable to that of 
PDA. Examination of the micrstructural data in these samples revealed significant 



63 


dilfcrcncc (as much as 20%) between plane dihedral angles in the two disector planes (Si 
and S 2 ) at same triple edges. This suggests that some of the triple edges must be curved 
(while the analysis developed here requires a straight edge) 

In Older to obtain the distribution of relative grain boundary energies from distribution of 
TDA, the following expression was used: 

r, _ r, 

sin a, sin a, sin«, 

where yi, yj and yi the energies of the grain boundary planes making a triple edge. Using 
equation 5.1, the ratio of the minimum energy to maximum energy: 


where yi> yj -’yi . 

figure 5.8 shows the distribution of the relative grain boundary energy ( — ) and the 

ri 

statistical parameter of the distribution are shown in table 5.1 (see figure 5.7). It could be 
seen that the grain boundary energy ratio increases towards unity with increasing 
annealing temperatures, which is expected, as grain boundary energies tend to equilibrate 
at higher temperatures. 
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Figure 5.5; TDA distribution of sample annealed at (a) 900°C for 60 minutes; 
(b) 900'’C for 180 minutes; (c) 950°C for 60 minutes;(d) 1000°C for 60 
minutes;(e) 1050°C for 60minutes 
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CHAPTER 6 

CONCLUSIONS 


The major conclusions of this study can be summarized as follows; 


1 . An automated methodology was evolved integrating image acquisition, 
processing and analysis to determine plane dihedral angles in polycrystalline 
materials. The algorithm developed uses morphological operations to clean 
microstructural images and performs thinning and pruning of grain 
boundaries. The automated procedure identifies triple points and the pixel 
coordinates of the corresponding grain boundary segments. The dihedral 
angles at each triple point are calculated by selecting the best-fit polynomials 
on the grain boundary segments on a sample of polycrystalline austenitic 
stainless steel. Parameters such as standard deviation of plane dihedral angles 
are compared as a function of annealing temperature. The results so obtained 
corroborated the efficacy of the scheme developed for evaluating dihedral 
angles. 

2. The standard deviations of plane dihedral angles decreases with increasing 
annealing temperatures, that is, the spread of the plane dihedral angle 
distribution reduces about the mean value of 120° class. This is expected 
because with increasing annealing temperatures grain boundaries tend to 
equilibrate. 

3. It is also evident that as the annealing temperatures were increased the 
standard deviations of true dihedral angles decreased and the peaks shifted 
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towards 1 20° . This explains that at higher temperatures the gram boundaries 
tend to equilibrate as also mentioned in the case of plane dihedral angles. 

4. The spread of the transformed true dihedral angles distribution is considerably 
less than the spread of the plane dihedral angles distribution. This has been 
attributed to the fact that the transformation from plane to true dihedral angles 
distribution results in the elimination of the variance due to random sampling 
of triple edges by the plane of polish. It has been shown that the parameters, 
such as, standard deviation and peak frequency of true dihedral angles 
distribution are much more sensitive to deformation as compared to the 
parameters for plane dihedral angles distribution. 

5. 'fhe standard deviations of the ratio of grain boundary energy decreased with 
increasing annealing temperatures and the mean value of the ratio 
progressively increased (values moved closer to 1 .0), as expected. 

6. The technique developed in this work to estimate the distribution of true 
dihedral angles from the measured PDA is not just limited to 316L austenitic 
stainless steel but can also be applied to any polycrystalline material. 

The suggestions and the scope for the future research are given below 


• To match the corresponding triple points to form a gram boundary edge 
the algorithm in appendix A.5 is used, the procedure was only partially 
successful and manual intervention was required for lower temperature 
annealed samples. This problem can be resolved by additionally 
comparing the orientation of corresponding grain boundary boundaries. 

. The study can be applied in the field of textures and sintered 


microstiuctures. 
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• To further increase the statistical accuracy of the disector technique the 
number of triple point analyzed should be increased to several hundred. 
This can be achieved by creating a montage consisting of perfectly 
matched adjacent images horn a plane section. 
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Appendix A 


Computer Programs 


A.1 Program to read image data and convert into coordinate 
system 

/* mod_exgb.c*/ 

^include <stdio.h> 

#include <conio.h> 

#define MAXPIXEL 15000 /* size for the max points analysed*/ 

#deflne M AXSEG 1 000 /* maximum pixels in a gb segment */ 

#define GBPIXEL 'O' /* value of gb pixel */ 

int getap(int value[][2], int nc, int n, int p[], int m); 
int test_tp(int value[][2], int p[]); 

main() 

{ 

int p[3],q[3],gb[MAXSEG],value[MAXPIXEL][2]; 
int i,j,k.,l,nc, count; 
char c; 

for(i=0,j=0,nc=0;(((c=getc(stdin))!=EOF) && (nc<MAXPIXEL));) 

{ /* reading pixel coordinates */ 

if(c=='\n') 

{ /* end of line */ 

/*cprintf("(%i,%i)\n",i,j);command*/ 
/*to print on the monitor*/ 

i=i++; 

j=0; 

} 

if(c == GBPIXEL) { 


value[nc][0] = j; /* array stores x,y position */ 

value[nc][l] = i; /* of gb pixel*/ 

nc++; 

} 

if((c == ’O’) II (c ==’!'))]++; 

} 

if (nc>-MAXPIXEL) cprintf(" Warning: Could not read all of image 
data!\n"); 

count = 0; /* set triple point counter */ 

for (i=0;i<nc;i++) 

{ 

k = getap(value,nc,i,p,3); 
if (k == 3) { /* k=3 : it is a triple point */ 

if(test_tp(value,p)==l ) continue; 
if (count > 0) printf("-l\n");/* end of line marker */ 
for (j=0y"=^3;j++) 

{ 

1 = 2 ; 

gb[0] = i; 
gb[lj = pO]; 
do { 

if ((k = getap(value,nc,gb[l-l],q,3)) = 2) { 
if(gb[l-2]==q[0]) /* exclude the duplicate pixel */ 

gb[l] = q[l]; 

else 

gb[l] = q[0]; 

1 ++; 

} 

} while ((k == 2) && (1 < MAXSEG)); 

for(k=0;k<l;k++) printf(”%i %i '’,value[gb[k]][0],value[gb[k]][l]) 
if (j < 2) printf("- l\n"); /* end of line marker */ 

} 

count++; 

} 

} 

printf(" %d\n", count); 
printf(”-99\n"); /* end of file marker */ 

} 



int getap(int value[][2], int nc, int n, int p[], int m) 

{ /* returns no. of adjacent pixels touching pixels n */ 

/* on return, p[] contains the first m adjacent points */ 
int i,j,k,x,y,dx,dy; 

X = value[n][0]; 
y = value[n][l]; 
i = n; 

j = -i; 

k = 0; 
do { 
i = i +j; 

if (i >= nc) break; /* no more gb pixels to check */ 

if(i<0) { 
i = n; 

j = i; 

} 

else { 

dy = value[i][l] - y; 
if((dy<-l)||(dy>l)){ 

if (j ==-!){ 

i = n; 

j = 1; 

} 

else 

break; /* no need to check further */ 

} 

else { 

dx = value[i][0] - x; 
if((dx>-2) && (dx<2)){ 

if (k < m) p[k] = i; /* found a pixel "i" touching pixel "n" */ 
k+4-; 

} 

} 

} 

} while(l); 
retum(k); 

} 



int testjp(mtvalue[][2],mtp0) 

int xi3],y[3]; 
int i,j; 

for(i=0;i<3;i++) 

^xfil=value[p[i]][0]; 

y[i]=value[p[i]][l]; 


for(i=0;i<3;i++) 

( 

j=i+l; 

if(j==3)j=0; 

( if(y[il==(yO]+^)llyW (yW- )) 

return! I ); 

} 

return! 1); 

} 

} 

return(O); 

} 



/K.2 Program for evaluating Plane dihedral angles 


/* modtpaSx*/ 

#include <errno.h> 

#include<stdio.h> 

#include<alloc.h> 

#include<math.h> 

#include<conio.h> 

#include "modcfitS .c” /* to include this file given in appendix A.3*/ 

#defme MAX_DATA 100 /* maximum number of data points per grain 

segment*/ 

#define MAX_TP 400 /*maximum number of triple points per image*/ 
#define DISTANCE(xl,yl,x2,y2) ((xl-x2)*(xl-x2)+(yl-y2)*(yl-y2)) 
typedef struct { 

int order, npts; /*order= order of polynomial fit; npts=no of pts on 

GB line */ 

float x,y,r2,*a; /* r^ = correlation coefficient*/ 

} polycoeff; 

typedef struct { 

float tx,ty,beta[3]; /* tx, ty = coordinates of triple points; beta= array 

of PDA*/ 

polycoeff g[3]; 

} triplepoint; 

int tpfit(triplepoint t[], int ntp); /* declaring five functions*/ 
int gb_line(float x[], float y[], int *ndat); 
void pda(triplepoint t[], int ntp); 

get_distance(double d[], double si, double s2, triplepoint t, int i); 
double slope(float a[], int order, int x); 


main() 

{ 


int ntp; 


/* 


} 


triplepoint t[MAX_TP]; 

cprintf("\nProgram Begin\n"); 

tpfit( "d :\\tc\\source\\dataout.txt" ,t,&ntp);*/ 

tpfit(t,&ntp); 

pda(t,ntp); 

cprintfC'ntp = %i\n",ntp); /* prints number of triple points 

evaluated*/ 


void pda(triplepoint t[],int ntp) 

{ 

int k,i,j,p,q,r; 
double s[3],d[3]; 

FILE *fp; 
char check; 


if ((fp = fopenC'coeff log";’w")) ==NULL) cprintf("unable to create 
coeff.logVn"); 

fprintf(fp,"#,tpx,tpy, order, a0,al,a2,a3,a4,a5,a6\n"); 


printf("#,tpx,tpy ,x_ 1 ,y_l ,x_2,y_2,x_3 ,y_3 ,ord 1 ,r2_l ,ord2,r2_2,ord3 ,r2_3 ,be 
tal ,beta2,beta3,sumbeta,x,y\n"); /*tpx, tpy coord of triple junctions*/ 

/*x_l,y_l= coordinates of point on first 
grain segment ; r2_l=correlation coeff of 
first grain segment*/ 

for(k=0;k<ntp;k++) 

{ 

if(t[k].g[0].npts<=4||t[k].g[l].npts<=4||t[k].g[2].npts<=4) 

{ 

/* printf("\n%f %f has too few points to continue\n",t[k].tx,t[k].ty);*/ 

continue; 

} 


if(t[k].tx=-l) 

continue; 



printf("%i,%f,%f,",k+ 1 ,t[k] .tx,t[k] .ty); 

for(i=0;i<3;i++) loop to calculate slope for each line*/ 

{ 

s[i]=slope(t[k].gfi].a,t[k].g[i].order,t[k].tx); 

/* prinlf(" slope %cl is %lf ",i,s[i]);*/ 

} 

/* printf("\n");*/ /* end of loop for slope*/ 


for(i 0;i<3;i++) /* loop to calculate distance from triple point 

to x_l and y_l */ 


{ 

j=i+l; 

ifG==3)j=0; 

get_distancc(d,s[i],sO],t[k],i); 

t[k].bcta[i|=acos((d[0]+d[l]-d[2])/(2*sqrt(d[0])*sqrt(d[l])))*(180/3.1412); 

} 

check = ' 

if((t[kl.beta[0]+t[k].beta[l]+t[k].beta[2])<360){ 

check = '*'; /* adding check marks to angles which are not 

adding upto 360 degrees in the first run*/ 

for(p=0;p<3;p++){ 
q=p+I; 
if(q==3) q=0; 

r=q+l ; 
if(r==3) r=0; 

if(((t[ k] .beta[p]+t[k] .beta[q])>=(t[k] .beta[r]- 
0.001))&&((t[k].beta[p]+t[k].beta[q])<=(t[k].beta[r]+0.001))) 

{ 

t[kj.beta[p]=l 80-t[k].beta[p]; 
t[k].beta[q]=l 80-t[k].beta[q]; 

} 

} 

for(p=0;p<3;p++) printf("%i,%f,",t[k].g[p].order,t[k].g[p].r2); 
for(p=0;p<3 ;p++) printf("%f,",t[k] .beta[p]); 
printf("%f',t[k].beta[0]+t[k].beta[l]+t[k].beta[2]); 



if (check == ' ') 
printf("\n"); 

else 

prinlf(",%c\n", check); 

for(p=0;p<3;p++) 

{ 

fpri ntf( fp;'%i,%f,%f,o/oi'',k+ 1 ,t[k] .tx^k] .ty,t[k] .g[p] .order)^ 
for(i-0;i<=t[k].g[p].order;i++)fprintf(^,",%f',t[k].g[p].a[i]); 

fprinlf(rp,"\n"); 

} 

} 

fclose(fp); 

} 

get_distancc(double d[], double si, double s2, triplepoint t, int i) 

{ 

double c 1 ,c2,y 1 ,y2,x 1 ,x2,ml ,m2,kl ,k2; 
intj; 

j=i+l; 

ifG==3)j=0; 

c I =t.ty-s 1 *t.tx; /* equation of straight line*/ 

c2=t.ty-s2*t.tx; 

xl=t.g[i].x; 

x2=t.gfjl.x; 

yl=t.g[i].y; 

y2=t.g[jj.y; 


if (si == 0) { 
xl =t.g[i].x; 
yl = t.ty; 

} 

else { 

ml=-(l/sl); 

kl=(yl-(ml*xl)); 

xl=((kl-cl)/(sl-ml)); 

yl=(((ml*cl)-(sl*kl))/(ml-sl)); 

} 

if(s2==0) { 



x2 = t.glj].x; 
y2 = t.ty; 

} 

else ( 

m2=-(l/s2); 

k2=(y2-(m2*x2)); 

x2=((k2-c2)/(s2-m2)); 

y2=(((m2*c2)-(s2*k2))/(m2-s2)); 

} 

d[0]=DISTANCE(t.tx,t.ty,xl,yl); 
d[ 1 ] =DI ST ANCE(t.tx,t.ty ,x2,y2) ; 
d[2]=DISTANCE(xl,yl,x2,y2); 
printfC'%lf,%lf,",xl,yl); 

} 


double slope( float a[], int order, int x) 

{ 

int i; 

double sum = 0; 


for(i=l ;i<=order;i++) 
sum=sum+i*((double) a[i])*pow(x,i-l); 

retum(sum); 


int tpfit(triplepoint t[], int *ntp) 

{ 

x[MAX_DATA],y[MAX_DATA],sig[MAX_DATA],a[MAX_SIZE] 

int ndat,i,ord,j,l,err,fe,upperord; 

float r2; 

/* FILE *fp;*/ 

/* DATA INPUT & CALCULATIONS */ 
for (i= 0 ;i<MAX_DATA;i++) sig[i] = 1; 



/* fp = (FILE *) fopen(fname,"r")- 
if(fp==NULL) { 

printfC'file error: %i\n",ermo); 
exit(O); 

} 

*/ 

(*ntp) = 0; /* counter for number of triple points*/ 

do 

{ 

for(i=0;i<3;i++) 

{ 

fe = gb_line(x,y,&ndat); 

err = polyfit(x,y,sig,ndat,l,6,&ord,a,&r2,l,x[0],y[0]); 
if (!err) 

{ 

t[*ntp].tx = x[0]; 
t[*ntp].ty = y[0]; 
t[*ntp].g[i]. order = ord; 
t[*ntp].g[i].x=x[ndat-l]; 
t[*ntp].g[i].y=y[ndat-l]; 
t[*ntp].g[i].r2 = r2; 
t[*ntp] .g[i] .npts=ndat; 
if ((t[*ntp].g[i].a = (float *) 
malloc((ord+l )*sizeof(fioat))) == NULL) 

{ 

cprintf("waming: out of memory in tpfit\n"); 
getchO; 

retum(-l); /* out of memory */ 

} 

forG=0;j<=ord;j++) t[*ntp].g[i].aD] = a^]; 

} 

else { 

t[*ntp].tx = -l; 
for (l=i+l;l<3;l++) 

fe = gbjine(x,y,&ndat); 

break; 

} 

} 

(*ntp)++; 

} while(fe != -99); 



return(O); 

} 

int gb Jine(floal x[], float y[],int *ndat) 

{ 

int num; 

(*ndat) = 0; 
do 
{ 

scanf( "%i " ,&num); 
if(num >= 0) { 

if(*ndat<MAX_DATA) 
x[*ndat] = (float) num; 
scanf("%i",&num); 
if (*ndat < MAX_DATA) { 
y[*ndat] = (float) num; 
(*ndat)++; 

} 

} 

} while(num >= 0); 
return(num); 

} 



A.3 Program for curve fitting 


/*modcfit3.c*/ 

#include <stdio.h> 

^include <stdlib.h> 

#include <math.h> 

#include cflt.h /*including header file cfit.h given in appendix A.4*/ 
/* 

main() 

{ 


float x[ 1 001, yf 1 00], sig[100], a[MAX_SIZE], xi; 
int i, ord; 

for (i=0;i<50;i++) 

{ 

xi = ((float) i)/50.0; 
x[il = xi; 

y[i] = 0.01+0.02*xi+0.04*xi*xi + ((float) random(100))/100; 
sigfil = 1; 

} 

polyfit(x,y,sig,50, 1 ,7,&ord,a,0,0,0); 


} 

*1 

int polyfit(float x[], float y[], float sig[],int ndat,int ordl,int ord2,int *ord, float 
a[], float *r2,int ianch, float anch_x, float anch_y) 

{ 

/* Polynomial Least Square Fit: 

X, y: "ndat" data points; sig: "ndat" standard deviations; 
ordl,ord2; search poly order ordl to ord2 and return best value in *ord (to 
implement) 

a: on return contains "*ord+r' polynomial coefficients 
ianch: non-zero value => set anchor at "anch_x,anch_j/’" 
return value: non-zero => no solution possible */ 



int temp, ij,dcgl,deg2, err, ia[MAX_SIZE]; 
float chisq,chisql,al[MAX_SIZE],ftest; 

double yc,ym; 

long double covar[MAX_SIZE][MAX_SIZE], a_x, st, sr; 

void (*fn)(float, long double [], int); /* function pointer */ 
if (ord2<ordl ) SWAP(ordl,ord2) 

if ((ordl<0)||(ord2>MAX_ORDER)) retum(-l); /* invalid polynomial 
order */ 

if (ndat<=ordl ) retum(-l); /* invalid no. of data points */ 
if (ndat<— ord2) ord2=ndat-l; /* reduce the upper polynomial order */ 

fn = polyfuncl; 

if (ianch) { /* set anchor point */ 

anchor_x = anch_x; 
anchor_jy = anch_y; 
fn = polyfunc2; 

} 

for( * ord==ord 1 ; * ord<=ord2 ;(* ord)++) 

{ 

for(i=0;i<=*ord;i++) ia[i] = 1; 
if (ianch) { /* set anchor for each order value */ 

ia[0] = 0; 
a[0] = anchor_ 3 ^; 

} 

err = gllsq(x,y,sig,ndat,a,ia,(*ord)+l,covar,&chisq,fii); 
if (err != 0) return(err);/* Error: solution not possible */ 
if (ianch) { 
a_x = 1 ; 

for (i=l ;i<=*ord;i++) 

{ 

a_x = a_x*((long double) anchor_x); 
a[0] -= a[ij*a_x; 

} 

} 

if((*ord)>ordl) 

{ 



degl=ndat-(*ord); 
deg2=ndat-( ( * ord)+ 1 ); 
if(chisql < chisq) break; 
if((dcg2!=0)&&(chisql !=chisq)) 

ftest-(((dcg 1 -deg2)*chisq)/((chisql-chisq)*deg2)); 

break; 
if(ftesl> 1) 
break; 

} 

for(i=0;i<=(*ord);i++) 
al|iha[ij; 
chisq l=chisq; 

} 

(*ord)=(*ord)-l; 

for(i=0;i<=(*ord);i++) a[i] = al[i]; 

ym = 0; /* calculation of correlation coefficient */ 

for(i=0;i<ndat;i++) ym += y[i]; 
ym = ym/ndat; 
sr = St = 0; 
for(i=0;i<ndat;i+4-) 

{ 

St += (long double) ((y[i] - ym)*(y[i] - ym)); 
yc = 0; 

for0=O;j<=*ord;j++) yc += a[j]*pow((double) x[i], (double) j); 
sr += (long double) ((y[i] - yc)*(y[i] - yc)); 

} 

if ((st == 0.0) II (sr == 0)) 

(*r2) = 1 ; 
else 

(*r2) = (float) ((st - sr)/st); 

/* printfC'tp %f %f Ord:%i Deg:%i Chi:%f r2: %f a: ",x[0],y[0],*ord,ndat- 
*ord-l,chisql,*r2); 

for(i=0;i<=*ord;i++) printf("%f ",al [i]); 
prmtf("\n"); 

*/ 

retum(O); 

} 



void polyfuncl (float x, long double afunc[], int ma) 

/* returns basis functions (atx) in "afunc" in apolynomial of order ”ma-l” */ 


int i; 

afunc[0] = 1 ; 
for (i=l ;i<ma;i++) 

{ 

afunc[ij = ((long double) x) * afunc[i-l]; 

} 

} 

void polyfunc2(float x, long double afunc[], int ma) 

{ 

/* returns basis functions (at x) in "afunc" in a polynomial of order "ma-1" 
under constraint that curve is anchored at: (anchor_x, anchor_y) */ 

int i; 

long double xl = 1, x2 = 1; 

afunc[0] = 1 ; 
for (i=l;i<ma;i+H-) 

{ 

xl = xl*x; 

x2 = x2*((long double) anchor_x); 
afunc[i] = xl - x2; 

} 


int gllsq(float x[J, float y[], float sig[], int ndat, float a[], int ia[], int ma, long 
double covar[l[MAX_SIZE], float *chisq, void (*gfunc)(float, long double 
[], int)) 

{ 

/* General Linear Least Square Routine: 

X, y: "ndat" data points; sig: "ndat" standard deviations; ma: # of 
parameters; 

ia: "ma" size with only zero or non-zero values; zero value implies, 
corresponding 

parameter held fixed to value stored in "a"; gfunc: user routine to supply 
"ma" 



basis functions evalated in array afunc (size- "ma")- 
returns ChiSquare in "chisq" and covariance matrix in covar 
return value: non-zero => no solution possible * / 

int err,i,j,k,l,ni,mfit=0; 
float ym,wt,sui'n,sig2i; 

long double beta[MAX_SIZE],afunc[MAX_SIZE]; 

for G=0;j<ma;j++) 
if (iaUl) mfit++; 

if (mfit == 0) return(-l);/* no parameters to be fitted */ 
for (j=0;j<mfit;j++) 

{ /* Initialize the (symmetric) matrix */ 

for (k=0;k<mfit;k++) covarO][k]=0.0; 
beta[jj=0.0; 

} 

for (i=0;i<ndat;i++) 

{ /* Loop over data to accumulate coeffcients of the normal equations 

(*gfunc)(x[i], afunc, ma); 

ym=y[i]; 

if (mfit < ma) 

{ /* Subtract off dependences on known pieces of the fitting function */ 

for (j=0;j<ma;j++) 
if (!iaO]) ym -= aO]*afuncO]; 

} 

sig2i=1.0/SQR(sig[i]); 
for (j=-l,l=0;l<ma;l+-+-) 

{ 

if (ia[l]) { 
wt=afunc[l]*sig2i; 
for (j++,k=-l,m=0;m<=l;m++) 
if (ia[m]) covarO][++k] += wt*afunc[m]; 
beta[j] += ym*wt; 

} 

} 

} 

for (j= 1 ;j<mfit;j++) /* Fill in above the diagonal from symmetry */ 
for (k=0;k<j;k++) 



covar[k] [j]=covar|j] [k]; 

if ((err=gaussj(covar, mfit, beta)) != 0) /• Matrix solution */ 

retum(err); /* no solution possible */ 
for (j=- 1 ,l=0;l<ma;l++) 

if (ia[l]) a[l]-beta[++j]; /* Partition solution to appropriate coeffecients a 

*/ 

*chisq— 0.0, /* and Evaluate ChiSquare of the fit */ 

for (i=0;i<ndat;i++) 

{ 

(*gfunc)(x[i ],afunc,ma); 
for (sum=0.0,j=0;j<ma;j++) 
sum += a[j]*afunc[j]; 

*chisq += SQR((y[i]-sum)/sig[i]); 

} 

covsrt(covar,ma,ia,mfit); /* Sort covariance matrix to true order of 
tting coecients */ 
return(O); 

} 

void covsrt(long double covar[][MAX_SIZE], int ma, int ia[], int mfit) 

{ 

/* Expand in storage the covariance matrix covar, so as to take into account 
parameters that are being held fixed (For latter, return zero covariances) */ 

int i,j,k; 

long double temp; 
for (i=mfit;i<ma;i++) 

for covar[i]0]=covar[j][i]=0.0; 

k=mfit; 

for(j=ma-l;j>=0;j--) { 

if(iaO]) { 

for (i=0;i<ma;i++) SWAP(covar[i][k],covar[i]0]) 
for (i= 0 ;i<ma;i++) SWAP(covar[k][i],covarD][i]) 
k-; 

} 

} 

} 

int gaussj(long double a[][]VLAX_SIZE], int n, long double b[]) 



{ 

/* Gauss-Jordon Method; uses Full Pivoting 
a: square matrix of size "n x n"; b: vector of size n 
on return: "a" contans inverse of "a" and the solution in "b" 
return value: = 0 => no error; 1= 0 => no solution possible V 

int indxc[MAX_SIZE],indxr[MAX_SIZE].ipiv[MAX SIZE]; 

int i,icol,irow,j,k,l,ll; 

long double big,dum,pivinv,temp,templ; 

for (j 0,j<n,j++) ipiv[j]— Oj /* set pivot tracking elements to zero */ 
for (i=0;i<n;i++) /* main loop for pivot search */ 

big=0.0; 

for (j=0;j<n;j++) 

if (ipivfj j != 1 ) /* skip, if diagonal element in jth row is a pivot */ 
for (k=0;k<n;k++) 

{ 

if(ipivfk|==0) 

{ 

if(a01[k]<O) 

tempi =-aO][k]; 

else 

tempi =aO][k]; 
if (tempi >=big) 

{ 

big = tempi; 
irow=j; 
icol=k; 

} 

} 

else if (ipiv[k] > 1) retum(l); /* Error: Singular Matrix-1 */ 

} 

} 

} 

++(ipiv[icol]); /* irow, icol should contain the position of pivot */ 



if (irow . icol) if irow icol, then pivot already in place */ 

r /I A 1 ^ \ cxtr f position pivot by swapping rows */ 

for (l=0;l<n;l++) SWAP(a[irow][l],a[icoliril) 
SWAP(b[irow],b[icol]) 


indxr[i j— irow, /* keep track of the original pivot position */ 

indxc[i]=icoI; 


if (a[icol][ico]] 0.0) retum(2); /* Error: Singular Matrix-2 */ 

pivinv— 1 .0/a[icol][icol]; /* Now divide the pivot row by the pivot 
element */ 

a[icol][icol]=l .0; /* ??? */ 

for (l=0;l<n;I++) a[icol][l] *= pivinv; 
b[icolJ *= pivinv; 

for (ll=0;ll<n;ll++) /* Reduce the rows except the pivot one */ 

{ 

if (11 != icol) 

I 

dum=a[lll{icol]; 

a[ll][icol]=0.0; 

for (l=0;l<n;l++) a[ll][l] -= a[icol][l]*dum; 
b[lll -= bficoipdum; 

} 

} 

} /* End of Main Loop */ 

for (l=n- 1 ;1>=0;1— ) /* Unscramble the solution by interchanging 

columns */ 

{ 

if(indxr[I] != indxc[l]) 
for (k=0;k<n;k++) 

SWAP(a[k][indxr[l]],a[k][indxc[l]]); 

} 

return(O); 

} 

n 



A.4 Header file called by curve fit program 

#define SWAP(a,b) {temp=(a);(a)=(b);(b)=temp;} 
static float sqrarg; 

#definc SQR(a) ((sqrarg=(a)) == 0.0 ? 0.0 : sqrarg* sqrarg) 

#defme MAX_ORDER9 

#definc MAX_SIZE MAX_ORDER+l 

static float anchor_x, anchor_y; 

int polyfit( float x[], float y[], float sig[],mt ndat,int ordl,int ord2,int *ord, float 

a[], float *r2,int ianch, float anch_x, float anch__y); 

void polyfunc 1 (float x, long double afuncQ, int ma); 

void polyfunc2(float x, long double afunc[], int ma); 

void covsrt(long double covar[][MAX_SIZE], int ma, int ia[], int mfit); 

int gllsq(float x[], float y[], float sig[], int ndat, float a[], int ia[], int ma, long 

double covar[][MAX_SIZE], float chisq[], void (*gfunc)(float, long double 

0, int)); 

int gaussj(long double a[][MAX_SIZE],int n, long double b[]); 



A. 5 Program to mako comparison of triple points on the two 
dissector planes 


/*compare.c*/ 

#include<stdio.h> 

#defme MAX 194 

typedef struct { long int x,y,dist; 

} typecoord; 
typecoord coord [MAX]; 
int sort( typecoord coord[],mt n); 
main() 

{ 


long int a[MAX][2],b[MAX][2],xl,x2,yl,y2; /* a and b array stores the 
values of tps from teh two sets */ 

FILE *fp; 
int i,j,k,n; 

fp=fopen("tpdata.txt","r"); 
i = 0; 

while(fscanf(fp," %ld %ld %ld %ld 
",&a[ilf01,&a[i][l],&b[i][0],&b[i][l])==4) 

{ 

printfC %ld %ld %ld %ld\n",a[i][0],a[i][l],b[i][0],b[i][l]); 
i++; 

} . 

n=i; 

for(i=0;i<n;i++) 

{ /* calculating the distance for each tp from every tp of the 

other set */ 


for(j=0;j<n;j++){ 
xl=a[i][0]; 
yl=a[i][l]; 
x2 = b[jl[0]; 
y2 = bG][l]; 

coord[j].dist = (((x2-xl)*(x2-xl))+((y2-yl)*(y2-yl))); 
coordOl-x = x2; 
coordOj.y = y2; 



} 

sort(coord,n); 

printfC %ld, %ld, ",xl,yl); 
for(k=0;k<5;k++) 

printf("%ld, %ld, ",coord[k].x,coord[k].y); 
printf("\n"); 

} 

fclosc(fp); 

} 


int sort(typecoord coord[],int n) 

{ 

long int i,j,tiTip,tmpl,tmp2; 
for (i=0; i<n-i; i++) { 
for G=0;j<n-l-i; j++) 

if (coord[j+l ].dist < coord[j].dist) { /* compare the two neighbors */ 
tmp = coordyj.dist; 
tmpl=coord[j].x; 

tmp2=coord[jJ.y; /* swap a[j] and a[j+l] */ 

coordyj.dist = coordy+lj.dist; 

coordyj.x=coordy+l].x; 

coordy].y=coordy+]].y; 

coordy+ 1 J.dist = tmp; 

coordy+l].x=tmpl; 

coordy + 1 j .y=tmp2; 

} 

} 

} 

L'l 



Appendix B 

SAMPLE INPUT/OUTPUT DATA OF DIGITAL IMAGES 
B.1 Sample bitmap of image 


mill 111 II on 111 in 1] 1111111111111111111011111111111111111111111111111111110111 

11 111] n 1 1 101 ^ I - 11 111 111 ] 111111111111111 11011111111111 iiiiiiiiiiii]_]_]_]_]_]_]_2iiioiii 
iiiimi]no]n]nin]iiiiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiiiii2_2^iiii]_iiiiiiiiioiii 
iiiiniinnoniiiiii] iiiiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiinioiii 
iiiiiininn ooooiiiiiiiiiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiinoiii 
1111111 nil 101111 000001111111111111111111011111111111111111111111111111111111001 
iiiniiinnom ni iiioooiiiiiiiiiiiiiiiioiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniiiio 
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iimiinioiniiiiiiiii ii iiiiiimiioimoiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiiniii 
iiiiiiiiioiini 111 imiiiiiiiimiiiioiiioiiiiiiiiiiimiiiiiiiiiiiiiiiiiiiiiiii 
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1110001 11 100111 11111111111111111111101111111110111111111111111111111111111111011 
000111 nn 11001 11111111111111111111011111111110111111111111111111111111111111011 
11111111 mi 13 001111111111111111110111111111111011111111111111111111111111111011 
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B.2 Triple Point Evaluation 


The above bitmap is subjected to the algorithm in appendix A.l, the triple points ( two 
triple points are bold highlighted) indcates x coordinate in column 1 and the y coordinate 
in column 2. The data is arranged with ascending order of the y coordinate. Each triple 
point end with a (-1 ) and the file terminates with a -99. The last number before -99 
indicate number of triple points evaluated from the image data shown in appendix B. 1 . 
The sample out put is given below; 


397 15 396 15-1 
397 15 398 15 -1 
397 15 397 16-1 
528 15 527 15 -1 
528 15 529 15 -1 
528 15 528 16 -1 
575 15 574 15 -1 
575 15 576 15 -1 
575 15 575 16-1 
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242 15 241 15 240 15 239 15 238 15 237 15 236 15 235 15 234 15 233 15 232 15 231 15 
230 15 229 15 228 15 227 15 226 15 225 15 224 15 223 15 222 16 221 17 220 18 219 19 
218 19 217 19 216 19 215 19 214 20 213 20 212 20 211 21 210 21 209 22 208 22 207 22 
206 22 205 22 204 22 203 23 202 23 201 24 200 24 199 24 198 24 197 24 196 24 195 24 
194 23 193 23 192 23 191 22 190 22 189 22 -1 

301 38 302 37 302 36 303 35 303 34 304 33 305 32 306 32 307 31 308 31 309 31 310 31 
311 31 31231 313 30 314 30-1 

301 38 302 39 302 40 303 41 303 42 304 43 305 44 306 45 307 46 307 47 307 48 307 49 
308 50 309 51 310 52 310 53 311 54 311 55 311 56 312 57 312 58 313 59 313 60 314 61 
314 62 315 63 315 64 316 65 317 66 318 67 319 68 320 69 -1 
520 43 519 43 -1 
520 43 520 42-1 

520 43 521 44 522 45 523 46 523 47 524 48 525 49 526 50 527 51 528 52 528 53 529 54 
529 55 530 56 531 57 532 58 533 58 534 59 535 60 536 61 537 62 538 63 539 64 -1 
401 47 401 46-1 
401 47 402 47-1 

401 47 400 48 399 48 398 48 397 49 396 49 395 50 394 50 393 50 392 51 391 52 390 53 

389 54 388 55 387 56 386 56 385 57 384 58 383 59 382 60 381 60 380 61 379 61 378 62 

377 63 376 64 375 65 374 66 373 67 372 68 371 69 370 69 369 70 368 71 367 72 366 73 

365 74 364 75 363 75 362 76 -1 

50 62 49 62 -1 
50 62 50 61 -1 

50 62 51 63 52 64 53 65 54 65 55 66 56 67 57 68 58 69 59 70 60 71 61 72 62 73 63 74 64 

75 65 76 65 77 66 78 66 79 67 80 67 81 68 82 69 82 70 83 71 83 72 84 73 85 74 86 75 87 

76 87 77 87 78 88 79 89 80 90 81 91 82 91 83 92 84 92 85 93 86 93 87 93 88 93 89 94 90 

94 91 94 92 94 93 94 94 94 95 95 96 95 97 95 98 95 99 96 100 96 101 96 102 97 103 98 

-1 

14 63 14 62-1 
14 63 15 63 -1 
14 63 14 64 -1 

539 64 538 63 537 62 536 61 535 60 534 59 533 58 532 58 531 57 530 56 529 55 529 54 
528 53 528 52 527 51 526 50 525 49 524 48 523 47 523 46 522 45 521 44 520 43 -1 
539 64 540 64 -1 
539 64 539 65-1 
325 66 325 65 -1 


325 66 326 66 -1 
325 66 324 67 323 67 
320 69 321 68 322 68 
320 69 319 68 318 67 
312 57311 56311 55 
306 45 305 44 304 43 

320 69 321 70 321 71 

321 81 321 82 321 83 
323 93 323 94 323 95 
323 104 323 105 323 

323 114 323 115 323 

324 124 324 125 324 
324 134 324 135 324 


322 68 321 68 320 69 -1 

323 67 324 67 325 66 -1 

317 66 316 65 315 64 315 63 314 
311 54 310 53 310 52 309 51 308 
303 42 303 41 302 40 302 39 301 
321 72 321 73 321 74 321 75 321 
321 84 322 85 322 86 322 87 322 
323 96 323 97 323 98 323 99 323 
106 323 107 323 108 323 109 323 
116 323 117 324 118 324 119 324 
126 323 127 323 128 323 129 323 
136 325 137 325 138 326 139-1 


62 314 61 313 60 
50 307 49 307 48 
38-1 

76 321 77 321 78 
88 322 89 323 90 
100 323 101 323 
110 323 111 323 
120 324 121 324 
130 323 131 323 


313 59 312 58 
307 47 307 46 

321 79 321 80 
323 91 323 92 
102 323 103 
112 323 113 
122 324 123 
132 324 133 
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B.3 Evaluation of PDA 

The data set of appendix B.2 is subjected to algorithm in appendix A.2. this out put is of 
section Si of the dissector, similar output of section Si is obtained. The parameters 
evaluated are highlighted below: 

# = triple point number 
tpXjtpy = triple point coordinate 

x_l,y_l,x_2,y_2,x_3,y_3 = arbitrary point on vectors r^,r 2 and respectively 
ordl,r2_l,ord2,r2_2,ord3,r2_3 = order and correlation coefficient of grain boundary 
segments 

betal,beta2,beta3 = three plane dihedral angles 
sumbeta = sum of plane dihedral angle at that triple point. 


#,tpx,tpy,x_l,y_l,x_2,y_2,x_3,y_3,ordl,r2_l,ord2,r2_2,ord3,r2_3,betal,beta2, betas, s 
umbeta,x,y 

4.189.000000. 22.000000.185.038903, 14.978055,288.200453,-6.520143,150.091494, 

78.284764.1.0. 958875.2.0.870281.5.0.998546.103.400429.140.712845.115.931717.360.0 
44991 

6.301.000000. 38.000000.195.128365, 56.698689,306.141395, 26.489728,319.239370, 

69.441249.4.0. 968997.3.0.977190.1.0.987728.124.100906.125.827988.110.116104.360.0 
44998 

13.320.000000. 69.000000.325.039516. 66.067918.292.341493.52.623289.354.090394, 

120.351985.1.0. 930579.3.0.993984.3.0.781635.119.193329.154.227707.86.623962.360.0 
44998 

14.141.000000. 72.000000.154.185542. 20.993644.174.423261.80.349695.108.831854, 

103.01 1990.3.0. 1 19484.5.0.998520.3.0.987638.89.543457.122.037247.148.464294.360 0 
44998 

16.221.000000. 82.000000.218.661283.78.650566.257.380944. 101.782187.216.267892, 

88.411804.2.0. 975546.2.0.993865.3.0.962158.153.478760.97.905502.108.660736.360.04 
4998 

22.545.000000. 105.000000.540.836580. 64.809744.5 17.500365, 

131.719387,587.473649, 

119.991877.2.0. 950238.2.0.991856.2.0.991766.128.277237.116.397568.115.370201.360. 

045006 


4 



300 . 713774 . 5 . 0 . 994832 . 2 . 0 . 994269 . 3 . 0 . 983531 . 119 . 494713 . 130 . 857193 . 109.693092 360 
044998 

52 . 455 . 000000 . 289 . 000000 . 457 . 056193 . 260 . 857976 . 346.016923 
330 . 016189 , 500 . 991958 , 

316 . 013692 . 2 . 0 . 713324 . 3 . 0 . 994060 . 1 . 0 . 984857 . 114 . 817245 . 128 . 963959 . 116 . 263786 . 360 . 
044991 

55 . 573 . 000000 . 299 . 000000 . 563 . 1 2743 1 . 286 . 89606 1 . 632. 1 99 1 3 6 , 
288 . 936322 , 559 . 848006 , 

338 . 292184 . 3 . 0 . 987681 . 2 . 0 . 991606 . 3 . 0 . 957694 . 119 . 569443 . 118 . 169228 . 122 . 306335 . 360 . 
045006 

57 . 78 . 000000 . 320 . 000000 . 53 . 404932 . 288 . 746559 . 44 . 177936 , 346 . 382090 , 59 . 670518 , 

337 . 303324 . 2 . 0 . 993361 . 4 . 0 . 771361 . 6 . 0 . 966311 . 90 . 234795 . 174 . 604034 . 95 . 161171 . 360.00 
0000 ,* 

58 . 343 . 000000 . 322 . 000000 . 249 . 349618 . 244 . 210223 . 457 . 150689 , 

300 . 328216 . 338 . 897731 . 367 . 721656 . 3 . 0 . 989624 . 4 . 0 . 993320 . 1 ,- 
2 . 0850 1 3 , 129 . 552002 , 1 05 . 890045 , 124 . 602959 , 360.045006 

61 . 25 . 000000 . 333 . 000000 . 42 . 832672 , 321 . 524222 , 46 . 136218 , 298 . 301244 , 32 . 788625 , 

355 . 025423 . 3 . 0 . 961198 . 5 . 0 . 843150 . 6 . 0 . 995228 . 154 . 106125 . 129 . 194504 . 76 . 699379 . 360.0 
00008 ,* 

64 . 107 . 000000 . 347 . 000000 . 10 . 882196 , 375 . 476277 , 137 . 577195 , 

271 . 559529 , 127 . 061884 , 

460 . 752726 . 2 . 0 . 991375 . 3 . 0 . 970328 . 2 . 0 . 985559 . 128 . 582230 . 147 . 952957 . 83 . 509819 . 360.0 
45006 

65 . 47 . 000000 . 350 . 000000 . 25 . 928435 . 351 . 081940 . 47 . 796537 , 328 . 028878 , 47 . 740974 , 

350 . 729626 . 5 . 0 . 996002 . 2 . 0 . 099828 . 3 . 0 . 943679 . 95 . 027489 . 132 . 498184 . 132 . 519333 . 360.0 
45006 

68 . 370 . 000000 . 410 . 000000 . 350 . 899564 , 365 . 205236 , 458 . 072684 , 
380 . 254347 , 361 . 370137 , 

422 . 128805 . 3 . 0 . 995838 . 1 . 0 . 964663 . 2 . 0 . 936394 . 94 . 443275 . 144 . 112427 . 121 . 489304 . 360.0 
45006 

70 . 219 . 000000 . 457 . 000000 . 170 . 61 1496 . 479 . 470269 . 278 . 635686 , 
337 . 799344 , 238 . 446494 , 

474 . 503956 . 3 . 0 . 968878 . 4 . 0 . 9921 15 . 1 . 0 . 969967 . 141 . 505173 . 105 . 425179 . 1 13 . 1 14647 . 360 . 
044998 

71 . 137 . 000000 . 459 . 000000 . 70 . 727878 , 375 . 941369 , 158 . 713392 , 

450 . 950264 , 123 . 726194 , 

485 . 806478 . 4 . 0 . 986472 . 4 . 0 . 984992 . 4 . 0 . 972689 . 108 . 258652 . 136 . 701477 . 115 . 084862 . 360 . 
044991 

73 . 448 . 000000 . 465 . 000000 . 409 . 984019 . 465 . 103027 . 470 . 013531 , 
462 . 180457 , 453 . 342319 , 

479 . 505634 , 3,0 967119 , 3 , 0 . 864450 , 1 , 0 . 627191 , 172 . 878052 , 77 . 090088 , 110 . 076866 , 360.0 
45006 


B.4 Comparison of Triple Points of Disector Sections 
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50 

265 

451 

222 

521 

159 

519 

151 

395 

134 

334 

156 

586 

278 

257 

235 

255 

245 

334 

156 

158 

288 

135 

260 

337 

335 

584 

252 

586 

278 

633 

300 

521 

159 

519 

151 

491 

370 

459 

261 

491 

370 

521 

159 

519 

151 

586 

278 

395 

134 

563 

287 

586 

278 

633 

300 

491 

370 

521 

159 

519 

151 

455 

289 

491 

370 

337 

335 

586 

278 

521 

159 

340 

380 

573 

299 

586 

278 

633 

300 

491 

370 

521 

159 

519 

151 

78 

320 

126 

317 

113 

274 

50 

265 

135 

260 

158 

288 

343 

322 

337 

335 

340 

380 

370 

427 

255 

245 

491 

370 

25 

333 

50 

265 

126 

317 

113 

274 

135 

260 

158 

288 

107 

347 

126 

317 

113 

274 

158 

288 

135 

260 

50 

265 

47 

350 

50 

265 

126 

317 

113 

274 

135 

260 

158 

288 

370 

410 

370 

427 

340 

380 

337 

335 

491 

370 

222 

459 

219 

457 

222 

459 

136 

458 

340 

380 

370 

427 

126 

317 

137 

459 

136 

458 

222 

459 

126 

317 

158 

288 

113 

274 

448 

465 

370 

427 

491 

370 

340 

380 

337 

335 

600 

600 
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